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Use  of  seismic  sensors  in  battle-field  surveillance  is  frequently  made 
difficult  by  burled  geophone  response  to  above  ground  acoustic  disturbances. 
ITie  response  to  such  acoustic-to-seismic  coupling  frequently  generates  signs 
in  the  same  frequency  range  and  of  the  same  amplitude  as  those  purely 
seismic  signals  that  are  of  interest  in  surveillance.  This  report  is 
concerned  with  the  theoretical  understanding  of  the  phenomenon  of  acoustic 
seismic  coupling. 
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1.2  In  the  work  reported  here,  the  penetration  of  sound  into  the  ground  is 
predicted  for  both  spherical  and  plane  acoustic  waves  incident  upon  two  models 
of  the  ground  viz  (i)  a  rigid  porous  solid  half-space  and  (il)  a  poro-elastic 
granular  half-space.  In  addition  the  plane  wave  analysis  is  extended  to 
incidence  on  a  rigidly-backed  finite  rigid-porous  layer 

1.3  In  the  frequency  range  of  interest  (1  -  1000  Hz)  the  rigid  porous  model  of  the 
ground  is  found  to  be  sufficiently  locally-reacting  to  predict  identical 
results  for  spherical  and  plane  wave  incidence.  The  predictions  are  found  to 
correspond  to  buried  microphone  measurements  in  sand  and  soils  of  above-ground 
acoustic  disturbances.  Furthermore  it  is  found  possible  to  predict  the  results 
of  model  measurements  using  continuous  and  pulsed  acoustic  sources  and 
microphone  receivers  above  and  below  the  surface  of  fibreglass  blanket. 

1.4  The  plane-wave  analysis  for  a  poro-elastic  half-space  predicts  the  existence 
in  general,  of  three  body  waves  i.e.  two  compressional  waves  which  we  call 
"fast"  and  "slow"  borrowing  from  the  literature  on  propagation  in  water 
saturated  sediments  and  one  shear  wave.  The  solid  particle  motion,  in  general 
is  the  result  of  the  contributions  from  all  three  wave  types.  At  normal 
incidence  the  contribution  is  mainly  from  the  fast-wave  type  travelling 
predominantly  in  the  solid  grains.  The  pore-fluid  particle  motion  has  a 
primary  contribution  from  the  slow  wave  type  and  for  dry  soils  is  practically 
identical  to  that  predicted  by  the  rigid  porous  model  of  the  ground. 

1.5  The  predictions  of  the  poro-elastic  models  for  normally- incident  plane-wave 
incidence  are  found  to  explain  qualitative  features  of  plots  of  intensity 
ratios  reported  in  the  literature  and  calculated  from  simultaneous, 
measurements  with  microphones  and  geophones  collocated  ie.  vertically  separated 
above  and  below  porous  ground  surfaces. 

1.6  The  analysis  for  spherical  wave  incidence  on  a  poro-elastic  half-space 
predicts  that  the  surface  motion  potentially  results  from  five  wave 
contributions  viz-  the  incident  wave  and  four  coupled  waves  in  the  poro- 
elastic  medium.  The  coupled  waves  consist  of  two  compressional  waves, 

a  shear  wave  and  a  Rayleigh-type  wave. 

1.7  Computations  of  the  near-  and  far-field  approximations  of  the  surface  motion, 
resulting  from  an  above-ground  acoustic  point  source,  for  a  horizontal  source- 
receiver  separation  of  100m  and  using  a  set  of  typical  soil  parameters,  show 
that,  for  source  heights  of  zero  and  1M  (grazing-incidence)  and  low-frequencies 
i.e.  less  than  100  Hz  the  Rayleigh-type  wave  is  dominant.  At  higher  frequencies 
for  grazing-incidence  or  for  all  frequencies  if  the  source  is  elevated  to  SOM, 
computations  reveal  that  the  direct-coupled  plane-wave  contribution  is  dominant. 
All  contributions  are  predicted  to  decrease  with  increasing  frequency  which  is 
in  accordance  with  the  trend  exhibited  by  maasured  data. 
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1  INTRODUCTION 

1.1  Scientific  interest  in  the  phenomenon  of  acoustic-to-seismic  coupling  has 
spanned  at  least  the  last  thirty-two  years  II.  1).  The  motivations  have 
included  seismic  exploration  using  sources  in  air,  studies  of  sonic  boom 
effects,  and  studies  of  the  ground- induced  motion  due  to  rocket  launches. 
Furthermore  ground-air  coupling  is  an  important  phenomenon  in  any  attempt 

to  monitor  sub-surface  conditions  acoustically  e.g.  to  detect  buried  objects 
using  sources  in  air. 

1.2  Systematic  experimental  studies  of  the  near  surface  seismic  signals  induced 
by  air-borne  acoustic  waves  are  relatively  recent.  Bass  et  al  [1.2,  1.3] 
have  used  continuous  signals  from  a  loudspeaker  source  and  impulsive  signals 
directed  along  a  line  array  of  triaxial  geophones  and  microphones  and  have 
reported  results  in  terms  of  the  squared  ratio  of  surface  seismic  particle 
velocity  to  sound  pressure  immediately  above  the  surface  termed  the  intensity 
ratio.  Similar  measurements  have  been  carried  out  by  Powell  [1.1]  and, 

recently, by  several  delegates  of  NATO  AC  Panel  III  RSG  II  on  mechanical  Haves  [1.4], 

1.3  Typical  results  show  that 

(i)  the  intensity  ratio  decreases  with  frequency  overall  but  has  several 
peaks  some  of  which  fall  outside  the  "exciting"  frequency  band 
(11)  the  coupling  is  not  very  sensitive  to  "incident"  angle  of  excitation 
(iil)  the  time  of  first  arrival  response  at  the  seismic  sensor  to  an  abcve-ground 
acoustic  impulse  corresponds  to  the  speed  of  sound  in  air  and  is  followed 
by  a  dispersive  Rayleigh  wave  train. 

(iv)  the  amplitude  of  seismic  disturbance  is  unaffected  by  depth  up  to  1.5  m. 

(v)  the  nature  of  the  acoustically- induced  seismic  disturbance  is  strongly 
affected  by  the  presence  of  snow. 


Cress  and  Flohr  (1.3)  have  attempted  an  analysis  of  the  coupling  using 
assumptions  of 

(i)  perfectly  flat,  homogeneous  and  infinitely  deep  ground. 

(ii)  continuous  vertical  displacement  across  the  boundary. 

(iii)  continuous  vertical  stress  across  the  boundary. 

(iv)  zero  horizontal  stress  components  at  the  boundary,  and 

(v)  compressional  wave  velocity  in  soil  -  34%  ms  1  (i.e.  and  equal  to 
that  in  air) . 

The  boundary  conditions  are  used  to  obtain  a  set  of  equations  relating  the 
reflected  and  transmitted  displacement  at  the  ground  surface  and  hence  the 
vertical  component  of  velocity.  It  is  assumed  further  that  the  acoustic 
pressure  can  be  related  to  the  particle  velocity  in  air  at  the  surface  by 
the  usual  plane  wave  relation. 

The  predicted  intensity  ratio  is  found  to  be  independent  of  frequency  and 
very  sensitive  to  angle  of  incidence.  Both  of  these  conclusions  are  at 
variance  with  measured  results. 

This  report  is  concerned  with  the  development  of  a  theory  for  the  phenomenon 
of  acoustic -to -seismic  coupling  which  is  in  accordance  with  observation.  Th 
starting  point  is  the  well-developed  theory  for  propagation  of  sound  from  a 
point  source  above  a  locally-reacting  boundary.  Appendices  1  and  2  repeat 
and  confirm  this  analysis.  In  Chapter  2  the  analysis  is  extended  to  predict 
the  field  within  the  locally-reacting  homogeneous  medium.  The  prediction  is 
found  to  depend  crucially  upon  the  acoustical  characteristics  of  the  locally 
reacting  medium.  Chapter  3  describes  "model"  experiments  carried  out  in  a 
anechoic  chamber  to  provide  data  for  detailed  comparison  with  theoretical 
predictions  and  Appendix  3  describes  the  basis  for  the  models  and  the 
instrumentation,  consequently,  in  Chapter  4,  a  development  of  the  Rayleigh 
model  (1.5)  for  such  a  medium  is  found,  in  combination  with  the  propagation 
(transmission)  theory,  to  give  predictions  that  agree  both  with  the  data 
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obtained  in  "model"  measurements  using  fibreglass  as  the  "ground" 
medium  and  with  microphone  probe  measurements  made  outdoors  in  sound 
and  soils. 


1.6  These  predictions  do  not  agree,  however,  with  the  geophone  data  reviewed 
in  1.3.  Thus,  in  Chapter  S  and  Appendix  4,  a  semi-infinite,  homogeneous 
poro-elastic  model  of  the  ground  is  adopted  (1.6)  and  used  to  give 
predictions  of  the  surface  intensity  ratio  after  a  review  of  model 
parameters  appropriate  to  sands  and  surficial  soils. 


Predictions  are  obtained  for  the  surface  intensity  ratio  due  to  acoustic 
plane  waves  normally  incident  an  both  a  semi-infinite  and  a  finite  depth 
rigidly-backed  poro-elastic  layer  and  obliquely-incident  on  a  semi -infinite 
poro-elastic  medium.  Finally,  since  the  plane  wave  analysis  will  not  be 
valid  near  grazing-incidence,  the  theory  of  spherical  wave  incidence  on 
an  alr/poro -elastic  boundary  is  derived  in  Chapter  6  and  Appendix  5  and  shows 
the  possibility  of  five  wave  contributions  at  the  Interface.  Near-  and  far- 
field  approximations  are  used  in  Appendix  5  to  show  the  relative  importance 
of  the  various  wave  types  in  acoustic -to -seismic  coupling.  The  predictions 
are  qualitatively  in  agreement  with  features  of  the  measured  data  that  are 
available  in  the  literature. 
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2  TRANSMISSION  OK  A  SPHERICAL  NAVE  AT  A  RIGID  POROUS  INTERFACE 
2.1  Introduction 

The  sound  field  due  to  a  point  source  above  a  plane  interface  separating 
two  semi-infinite  fluid  media  is  a  problem  of  great  interest  to  those 
concerned  with  the  propagation  of  sound  over  different  types  of  ground 
surface,  where  the  soil  can  be  approximated  by  a  rigid,  porous  absorbing 
medium  [2.1].  Considerable  effort  has  been  made  to  evaluate  the  field 
in  the  upper  medium  (air)  in  the  short-wave,  asymptotic  limit.  A  geometric 
approximation  is  adequate  for  sufficiently  large  elevation  of  the  source 
and/or  receiver,  but  this  breaks  down  at  grazing  incidence  where  small 

correcting  terms  play  an  increasingly  dominant  role.  The  grazing- 
incidence  case,  of  course,  is  precisely  that  which  applies  in  long-range 
propagation  studies,  (hie  of  these  correcting  terms  has  the  form  of  a 
surface  wave  and  in  Appendix  2  we  show  that  it  corresponds  to  a  physically- 

real  surface  wave. 

The  evaluation  of  the  field  in  the  lower  medium  has  recieved  comparatively 
little  attention,  but  Paul  ((2.2],  (2. 3])  and  Brekhovskikh  [2.4* ] 
discuss  the  problem  In  some  detail.  Our  motivation  for  considering  this 
case  arises  from  a  desire  to  compare  the  magnitude  of  sound  waves  penetrating 
into  the  ground  with  that  of  any  accompanying  seismic  disturbance. 

We  start  in  a  similar  manner  to  Paul  and  Brekhovskikh,  that  is  by  expressing 
the  solution  as  an  inverse  Hankel  transform.  The  resulting  integral 
is  then  evaluated  asymptotically  using  two  steepest  descent  contours  which 
pass  around  the  two  branch  points  of  the  integrand.  The  presence  of  a  pole, 
which  affects  the  solution  at  grazing  Incidence  over  ground  having  a 
high  impedance,  ie  also  taken  into  account.  Wa  only  consider  cases  for  which 
the  refractive  index  has  a  magnitude  greater  than  unity. 
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In  his  first  paper  [2.2]  Paul  integrates  around  two  vertical  (in  the  | 

plane  of  complex  wavenumber)  paths  originating  at  the  branch  points, 

as  opposed  to  the  true  saddle  paths  that  we  use.  This  initial  simplification  j| 

still  results  in  complicated  algebraic  espresslons  and  furthermore  no 

ready  physical  interpretation  is  possible  in  terms  of  geometric  rays.  j 

Be  deals  with  the  pole  by  subtraction,  and  while  this  method  appears  to 

be  essentially  equivalent  to  ours,  we  feel  that  the  final  result  is  not  1 

expressed  so  siaply.  | 

Paul  later  Considers  alternative  Integration  paths  [2.3]  which  in  | 

part  resemble  more  closely  those  used  here.  Be  does  not  consider  the 
presence  of  the  pole  in  this  work. 


Brekhovskikh  [2.4]  does  not  describe  his  derivation  in  detail,  but 
it  would  appear  that  he  uses  the  true  saddle  paths  as  he  is  able  to 
identify  geometrically  the  two  wave  types.  It  is  not  made  clear, 
however,  how  to  deal  with  a  complex  refractive  index  (i.e.  an  absorbing 
medium)  except  in  an  approximate  form,  and  the  problem  of  the  pole  is 
only  mentioned  briefly. 

In  a  more  recent  paper  by  Candel  and  Crance  [2.5]  the  related  problem  of 
a  horizontal  line  source  is  discussed.  Their  stationary  phase  method 
can  only  give  a  single  term  which  is  the  same  as  that  obtained  using 
geometrical  acoustics,  and  is  therefore  not  suitable  at  grazing  incidence. 
They  do  not  acknowledge  the  presence  of  a  lateral  wave,  nor  is  the  pole 
taken  into  account. 
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2.2  The  integral  representation  for  the  field 


Expressed  as  an  inverse  Henkel  transform  the  velocity  potential  in 
the  lower  medium  is  given  by  [2.2] 


♦ 


VUj.z^r) 


-Vizi  -  VjZ2  +  iKr 
ft  _ _ 

vx  +  MU2 


1 

H„(Kr)e 


iKr 


Kdx 


(2.1) 


where  z} ,  z^  are  the  depth  of  source  and  height  of  receiver  and  r  is 
their  horizontal  separation, 


Vi 


(K1 


k1)' 

l 


v*  -  (K  -  kj)*,  Rev  Rev2  >  0, 


k _,  k?are  the  wavenumber  (complex  in  general)  in  the  upper  and  lower 
media,  and 


N  «  Pi/P2  ,  the  ratio  of  densities  (again  complex) 


Equation  (2.1)  can  be  used  to  obtain  an  expression  for  the  field 
in  the  upper  medium: 


eikiRi 

4tRj 


#ikiRi 

4«R2 


(H,0,r) 


(2.2) 


where  B  is  the  sum  of  source  and  receiver  heights  and  Rj  ,  Rj  are  the 
distances  from  the  receiver  to  the  source  and  to  its  image. 
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The  steepest  descent  method  with  a  Dole  near  a  saddle  point 

The  rapid  oscillation  of  the  integrand  in  equation  (2.1)  for  large 
source-receiver  separations  can  be  handled  satisfactorily  by  using 
the  saddle-point  method. 

In  general  we  need  to  evaluate  integrals  of  the  form 


(2.3) 


in  which  the  behaviour  of  f (K)  causes  the  integrand  to  oscillate  rapidly. 
The  saddle  points  occur  where  V  (K)  *  0 ,  and  the  steepest  descent  paths 
emanate  from  these  in  a  direction  in  the  complex  plane  in  which  the 
exponential  term  decays  most  rapidly.  The  imaginary  part  of  f(K)  is 
constant  along  a  steepest  descent  path  and  can  be  taken  outside  the 
integral,  leaving  a  purely  real  exponential  term  which  decays  on  either 
side  of  the  saddle  point.  If  f 1 (K^)  «  0,  then 


f (K)  ■  f (K  )  -  S'  for  real  S  gives  the  behaviour  of  the 
8 

♦ 

exponent  along  the  steepest  descent  path.  The  integral  along  this  path 
becomes 


dX  -S’ 

A(K)  —  e  dS 


(2.4) 
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If  A(K)  —  can  be  expanded  in  a  power  series  in  s  with  a  large  radius 
cfo 

of  convergence  we  can  integrate  the  successive  terms.  With 
A(k)dK/dS  -  MS),  this  results  in 


SOP 


,ief(K*) 


[♦(0, 


(0)  + 


■] 


(2.5) 


If  one  or  more  (but  not  necessarily  all)  of  the  steepest  descent  paths 
are  added  in  such  a  way  that  the  resulting  contour  can  be  deformed  into 
the  real  axis  without  crossing  any  branch  points,  then  the  resulting 
integral,  along  with  the  contributions  from  any  poles  crossed,  will  be 
equal  to  the  integral  in  equation  (2.3) 


If  the  radius  of  convergence  of  the  series  for  Ms)  is  limited  by  the 
presence  of  a  pole  equation  (2.5)  may  be  inapplicable.  In  this  case 
we  write 


«(S)  -  <MS)  /  (S  -  sp> 


(2.6) 


and  use  the  series  for  MS/,  which  has  no  pole,  instead.  We  now 
obtain 


I  -  .*«•>  *  (o)I- 

SDP  L 

+  *’(0>  (Spl,  ♦  /?  ) 

+  *♦-  (0)  (8  *1.  +  S/*  ) 

P  P 

♦  1/6  W"'<0>  IS  *1.  ♦  S  */»♦  »/iT  ) 
P  P 
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where 


-S’ 
e _ 

s-s_ 


dS 


Noting  that  (MSp)  «  *><0)  +  *’(0)Sp  +  Ji|i"<0)SpJ  + 
and  making  use  of  equation  (2.6)  gives 


»  n 

SDP 


Jef(K»)  t*(Sp)  <I„  ^  ^  +  '2g~T  ♦  ••••> 


+  *(0)  +  j  *"(0)  + . ) 


(2.8) 


It  can  be  seen  that 
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I 

I 

I 
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] 


I  +  sief(KB>  (*'(0)  ♦  1/12*"M0)  ♦ . ) 

as  S  0. 

P 

The  expression  therefore  remains  finite  as  the  pole  approaches  the  saddle 
point. 

He  can  write 

I„  «  is  {*<Sp>  -  2e-SP  B  l-Im(Sp)}} 

where  W (S  )  »  e'V  erfc  (-iS  )  and  B[  ]  is  the  Beaviside  step  function, 
p  P 

This  second  term  is  cancelled  by  the  contribution  from  the  pole  residue 
if  the  pole  lies  between  the  SDP  and  the  real  axis,  leaving  the  total 
contribution 

I*  -  ^*(K*)  t*<Sp)  U*F  wCSp)  ♦  ^  +  •*‘) 

♦  4(0)  +  ^"(0)  ♦ - 1 


W 

* 

j 

i 


i 


(2.9) 
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To  Investigate  the  behaviour  of  this  expression  for  large  values 
of  3^  we  use  the  asymptotic  expansion 


I*  -  S*  ef(Ka)  [♦(0)  +  ^"(0) 

♦  2B  [-m(Sp)]  e-SP%  (Sp)l  (2.10) 


As  expected,  this  expression  approaches  the  behaviour  of  equation  (2.5) 
along  with  the  contribution  from  the  pole  when  it  lies  between  the  SDP 
and  the  real  axis.  In  this  application  the  latter  term  gives  rise 
to  the  faeous  'surface  wavs'. 


Thus  the  SDP  Integral  eay  calculated  using  equation  (2.5)  and  if 
the  pole  is  near  the  saddle  point  we  simply  add  the  extra  terms  shown 
in  equation  (2.9).  This  straightforward  procedure  appears  to  be 
much  less  cumbersome  than  previous  methods  such  as  subtraction  of  the 
pole. 

The  subsequent  means  used  to  evaluate  the  expression  in  equation 
(2.5)  are  well-established  and  Brekhovskikh  (2.4)  gives  a  good 
account  of  these,  and  the  appropriate  formulae. 


For  the  extra  tern  in  equation  (2.3.9)  we  use 


*<S  )  -  lim  (K  -  K  )A(K) , 


K-*K 


P 


(2 


where  is  the  position  of  the  pole  in  the  complex  K  plane. 


Returning  to  equation  (2.1)  we  have 


f(K)  »  -Vx*!  -  v*z2  +  iKr  (2 

If  —  ilfr 

and  A(K)  -  [H;  (Kr)e  ]  K/(Vi  +  Kv2) .  (2 

For  grazing  incidence,  with  Zj,  z 2  «  r,  there  are  two  saddle  points,  one 
near  K  *  k^  and  another  near  *  «  k^ .  The  steepest  descent  paths  through 
these  points  can  be  used  to  close  the  contour  of  integration  without 
enclosing  the  branch  points.  The  pole,  given  by  v2  +  Mv2  •  0,  can 
closely  approach  the  first  of  these  saddle  points  at  grazing  incidence 
if  M  1 ,  kj  ^  k j ,  and  will  be  dealt  with  accordingly  using  equation 
(2.9) 

At  sufficiently  large  source  heights  or  receiver  depths  closure  of  the 
path  of  integration  can  be  effected  by  using  only  one  of  the  two  paths. 

Ne  will  not  pursue  this  natter  further,  but  both  Paul  (2.3]  and 
Brekhovskikh  (2.4]  consider  such  cases. 
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The  geometrical  wave 


The  contribution  from  the  saddle  point  near  the  K  ■  k^  branch  point 
corresponds  to  the  ray  refracted  according  to  Snell's  law  in  the 
geometrical  approximation. 


We  locate  this  saddle  point  using  the  in ter at ion 


®i+i  ”  /<r  "  *2®i  /Yi> 


where  *  -ikia,  V2  *  -ikiy*  K  *  kx0  with  ft.  *  (1  -  a^)  , 


(o^  ♦  n*  -  1)*,  n  «  ki/kx. 


For  large  angles  of  incidence  (i.e.  *x  <  *)  Oi«  0  is  a  suitable 
starting  value,  and  the  calculations  may  be  performed  for  complex 
kj  and  k2  .  If  z,  *  0  then  o  -  0,  but  equation  <2.5)  not, 

becomes 


♦  yr*  +  8r) 


f  ilY  L,  *  1  fIV  .  1_  (f"')J 

^f-/  |*u  8  WT*  24  (r'F 

*  TrF  1  f  J 


(2.15) 
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2.5  The  lateral  wave 

To  find  the  second  saddle  point  we  take  into  account  the 
fact  that  at  grazing  incidence  it  is  close  to  the  K  =  k, 
branch  point,  and  consequently  |v2jis  small,  and  vx  is 
predominantly  real.  We  re-define  a,  6,  Y  and  use  the 
iteration 

Yi+  ■*  zi8i/(r  *  iziBj/aJ  (2.17) 

where  Vi  “  k,a,  V,  =  -ik,  Y»  K  •=  k,  B 
and  “  (1  “  Y^1  )  ^  / 

a,  -  a  -  -  yj>*.  r0  -  o 

This  gives 

f(K8)  »  -  kjOZj  +  i*,  (Yz»  +  Br)  • 

From  this  we  can  see  that  the  exponential  term  decays  as 

the  height  of  the  source,  increases.  Further  exponential  decay 

occurs  if  k,  has  a  significant  imaginary  part. 

The  steepest  descent  method  may  now  be  applied  in  the  usual 
way,  and  in  this  case  it  will  not  be  necessary  to  take  into 
account  the  presence  of  the  pole. 

Special  cases 

(i)  The  geometric  approximation. 

If  n,  the  refractive  index,  is  real,  then  equation 
(2.14)  is  satisfied  if 
a  ■  cos  0^, 

3  -  sin  9i, 
y  «  n  cos  ©r. 
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To  calculate  the  further  terns  In  equation  (2.9)  we  need  to 
locate  the  position  of  the  pole  using  vj  +  K\>2  =  0.  At  this 
point 
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sin  ©jVsinQ^  -  n,  (Snell's  Law) 
and  22 •  r  are  as  depleted! 


Taking  the  first  tern  in  equation  2.3.15  and  using  the 
first  tent  of  the  expansion  for  the  Hankel  function  so 


I-kxSt-f;,*  g 

x  expiki  (out!  *  Y**  +  Br) 

2  /sinSi  exp  iki  (*i/cos0i  +  n*2/cos0 

Ir  (iji/cos *0^^  ♦  «*/n  cos,0r)}^(n  cosB^  ♦  co«£) 


(2.18) 


which  is  the  geometric  approximation  given  by  Brekhovskikh 
ii)  Grazing  incidence  with  a  large /complex  refractive  index. 


If  n  has  an  appreciable  imaginary  component  then  the 
lateral  wave  will  decay  to  negligible  proportions  at 
large  distances. 


If  a,  is  sufficiently 
zaS/Y  «  r 


all  so  that  in  equation  (2.3.14) 


(2.19) 


Then  o  ■  *i 

r6 


i 
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Since  |8|  <  1/  a  is  small  at  grazing  incidence 
equation  (2.3.29)  becomes 

Z2 

-  «  1 

r(n2  -1)T 

Under  these  circumstances  01  «  cosO  where  tanO  * 
r/zj 


(2.20) 


If  in  addition  cosG  is  sufficiently  small  so  that  the 
terms  containing  Z2  are  negligible  in  f11 ,  f111,  flV 
Then  v(*l .  «r)  -*• 

V(2l#  o,  r) 

=  eiklZ2(n"  -  8in20  ’*  vUx,  0,  r) 

which  gives  a  simple  exponential  relationship  with  the  field 
on  the  surface  above  the  receiver.  This  expression  is 
derived  somewhat  less  rigorously  by  Brekhovskikh. 
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POINT  TO  POINT  PROPAGATION  MEASURMENTS  USING  ACOUSTIC  MODELLING 


Introduction 

In  order  to  validate  the  foregoing  theory  for  both  reflection  and 
transmission  of  spherical  waves,  measurements  of  point  propagation 
over  plane  absorbing  and  rigid  boundaries  were  made  in  the  small 
anec*»oic  chamber.  Two  cases  have  been  examined  in  detail: 


(i)  attenuation  over  a  rigid  boundary; 

(ii)  attenuation  over  and  within  a  porous  boundary. 

The  two  materials  selected  to  model  these  surfaces  were  varnished 
blockboard  and  Crown  150  fibreglass  quilt.  The  basis  and  facilities  used  for 
these  measurements  are  outlined  in  Appendix  3. 

Rigid  surface  measurements 

The  rigid  boundary  measurements  were  made  in  order  to  examine  the 
effects  of: 

(i)  source  height  uncertainty; 

(ii)  spectral  smoothing; 

(iii)  near  surface  turbulence;  and 

(iv)  source  and  receiver  size. 


The  source  and  receiver  were  placed  100  cm  apart  above  the 
varnished  blockboard  surface  as  shown  in  Figure  (3.1).  The  source 
and  receiver  heights  could  be  varied  independently. 

An  air  jet  was  used  as  the  sound  source  and  a  1/4"  microphone  as 
the  receiver.  A  block  diagram  of  the  measurement  system  is  shown 
in  Figure  (3.2) 

Source  height  errors 

For  a  rigid  boundary  the  excess  attenuation  minima  result  from  the 
destructive  interference  of  the  direct  and  reflected  waves.  The 
frequencies  at  which  these  minima  occur  are  therefore  very  sensitive 
to  source-receiver  geometry. 

Figure  (3.3)  shows  a  comparison  of  the  measured  and  predicted  excess 
attenuation  for  propagation  over  a  rigid  boundary.  The  prediction 
is  based  upon  a  source  and  receiver  height  of  5  cm.  It  can  be 
seen  that  the  location  of  the  predicted  attenuation  minima  do  not 
coincide  with  the  measured. 

The  source  height  estimate  was  based  upon  a  measurement  of  the 
distance  between  the  rigid  surface  and  the  top  of  the  air  jet 
nozzle.  However,  In  view  of  the  large  mass  of  air  flowing  through 
the  jet,  the  acoustic  centre  of  the  source  is  unlikely  to  coincide 
with  the  top  of  the  nozzle. 
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A  better  estimate  of  the  source  height  can  be  obtained  by 
varying  the  height  used  in  the  pure-tone  predictive  algorithm 
until  the  locations  of  the  measured  and  predicted  minima 
coincide.  The  variation  in  the  minima  location  and  source  height 
can  be  seen  in  Figure  (3.4)  increasing  the  source  height  increases 
the  path  length  difference  and  hence  moves  the  minima  to  lower 
frequencies. 

The  results  of  this  procedure  indicated  that  the  source  height 
is  actually  5.3  cm  above  the  rigid  surface.  A  comparison  of 
the  measured  and  predicted  excess  attenuation  curves,  based 
upon  a  source  height  of  5.3  cm  is  shown  in  Figure (3.5) 

.1.2.  Spectral  smoothing 

The  predicted  excess  attenuation  curves  shown  in  Figure  (3.3) 
and  (3.5)  exhibit  much  deeper  minima  than  the  measured  results. 

Part  of  this  discrepancy  arises  from  the  fact  that  the  predictions 
are  based  upon  pure  tone  interference,  whereas  the  measured  values 
are  based  upon  the  average  energy  in  a  250  Hz  wide  filter.  The 
effect  of  increasing  the  analysis  bandwidth  is  to  increase  the 
sound  pressure  at  the  minima  and  so  smooth  the  excess  attenuation 
spectrum. 

Correction  procedures  for  the  pure-tone  calculations  have  been 
derived  [3.1]  but  they  cannot  account  for  the  observed  differences. 
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Near  surface  turbulence 

Another  factor  to  consider  when  trying  to  account  for  the  observed 
difference  between  the  measured  and  predicted  excess  attenuation 
curves  is  the  effect  of  near  surface  turbulence. 

Daigle  [3.2]  has  proposed  modelling  turbulent  air  motion  as  an 
additional  propagation  path  length,  with  random  direction  and 
attenuation.  The  extra  path  length  disturbs  the  strong  interference 
pattern  due  to  the  surface.  A  correction  procedure  based  upon 
the  measurement  of  four  meterological  parameters  shows  good 
agreement  with  experimental  results.  These  parameters  are: 

-  variance  of  sound  velocity  deviations 

-  variance  of  air  temperature  deviations 
Ly  -  sound  velocity  correlation  length 

Lt  -  temperature  correlation  length 

Ly  and  l_t  are  measures  of  the  autocorrelation  function  of  the  wind 
and  temperature  variations.  The  time  delay  corresponding  to  maximum 
correlation  is  converted  to  an  equivalent  path  length  difference 
by  multiplying  the  velocity  of  sound  by  the  time  delay. 

Whilst  it  has  not  been  possible  to  measure  the  parameters  required 
to  correct  for  turbulence  effects,  an  examination  of  Daigle's 
results  indicate  that  they  could  account  for  the  observed 
discrepancies  between  the  measured  and  predicted  excess  attenuation 
curves. 
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Fibreglass  Measurements 


Two  sets  of  measurements  were  made  using  a  fibreglass  quilt  material  to 
model  a  porous  half-space,  the  first  to  measure  the  amplitude  ratio  of 
the  direct  and  transmitted  sound  fields  and  the  second  to  measure  the 
phase  difference  between  the  arrivals  of  those  fields 

Amplitude  Ratio  Measurements 

The  source  receiver  geometry  for  these  measurements  is  shown  in  Figure  (3.1) 
The  surface  receiver  was  a  standard  1/4"  microphone  with  a  normal  protection 
grid,  whilst  for  the  sub-surface  receiver  the  1/4"  microphone  was  fitted  with 
a  nose-cone.  The  nose-cone  serves  two  purposes,  first  its  shape  enables  one 
to  ease  the  microphone  through  the  fibreglass  quilt  and  second  it  improves 
the  directivity  response  of  the  microphone. 

The  sound  source  used  for  these  measurements  was  a  mid-range  loudspeaker, 
with  the  cross-over  network  removed,  and  the  source  signal  was  derived  from 
a  swept  frequency  oscillator  covering  the  range  100Hz  to  15kHz. 

In  order  to  minimise  the  effects  of  phase  distortion  that  can  arise  from 
incoherent  averaging  of  swept  frequency  signals,  a  time  gating  system  was 
employed  to  synchronise  the  source  signal  sweep  and  the  digital  sampling  of 
the  received  signals. 

The  narrow  band  spectrum  at  each  of  the  receiving  microphones  was  measured 
with  the  Nicolet  6606  spectrum  analyser,  and  the  amplitudes  ratio  obtained 
from  the  modulus  of  the  'transfer  function'.  Three  sets  of  results  are 
reproduced  here  as  Figure  (3.6)  Part  (a)  shows  the  amplitude  ratio  and 
phase  difference  for  a  source  height  of  5cm  and  a  receiver  depth  of  3cm. 
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The  ripple  in  the  transfer  function  at  8kHz,  arises  from  resonance  of 
the  loudspeaker  diaphragm.  Parts  (b)  and  (c)  show  the  results  for 
receiver  depths  of  7cm  and  10cm  respectively. 

A  test  of  the  reliability  of  these  measurements  was  performed  by  calculating 
the  coherence  function.  The  results  for  the  three  receiver  depths,  shown  in 
Figure  (3.7) .indicate  that  only  the  estimates  of  the  amplitude  ratio  at  a 
depth  of  10cm  for  frequencies  above  8kHz  may  be  erroneous;  since  the  value 
of  the  coherence  function  drops  below  0.8.  The  reduction  in  coherence  arises 
from  the  high  level  of  attenuation  at  this  receiver  depth. 

All  these  results  show  that  there  is  little  attenuation  of  the  transmitted 
sound  at  low  frequencies.  At  high  frequencies  the  attenuation  is  dependent 
upon  the  depth  of  the  buried  receiver.  At  10kHz,  the  upper  limit  of  the 
measurements,  the  attenuation  increases  from  16dB  at  a  depth  of  3cm,  to 
approximately  50dB  at  a  depth  of  10cm. 


.2.2.  Pulse  measurements 

The  time  delay  measurements  were  made  using  the  same  source  receiver 
geometry  described  in  section  (3.2.1 )%  but  with  the  loudspeaker  replaced 
by  an  electrical  spark  discharge  source.  This  device  produces  a  low 
intensity  sound  pulse  of  approximately  200  microseconds  duration. 
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The  trigger  transformer  circuit  used  to  generate  the  high  voltage  spark 
can  cause  significant  errors  due  to  electro-magnetic  radiation,  which 
swamps  the  response  of  the  condenser  microphone  capsules.  In  view  of 
this  the  spark  was  operated  at  its  minimum  energy  level  of  4mJ,  and  time 
averaging  used  to  improve  the  signal  to  noise  ratio. 

Examples  of  the  direct  (surface  microphone)  and  transmitted  pulses  are 
shown  in  Figure  (3*8)  for  a  source  height  of  5cm  and  a  receiver  depth 
of  3cm.  The  total  duration  of  the  signal  samples  is  4  milliseconds. 

The  direct  microphone  signal  shows  a  distinct  pulse  at  1  millisecond, 
corresponding  to  the  spark  discharge.  The  other  fluctuations  at  2.5 
and  3.5  milliseconds  correspond  to  the  arrival  of  reflections  from 
metallic  fittings  within  the  anechoic  chamber.  The  transmitted  pulse 
occurs  at  the  same  time,  but  the  trace  shows  no  subsequent  arrivals  due 
to  extraneous  reflections  within  the  chamber. 

The  spectrum  of  the  direct  and  transmitted  pulses  are  shown  in  Figure 
(3.9)  for  the  range  125Hz  -  50kHz.  The  upper  trace  corresponds  to  the 
direct  pulse  and  the  lower  trace  the  transmitted  pulse.  The  raggedness 
of  these  curves  is  due  partly  to  the  presence  of  reflection  and  partly 
to  the  electrical  noise  in  the  microphone  pre-amplifiers. 
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The  amplitude  ratio  and  phase  difference  for  the  direct  and  transmitted 
sound  fields  are  shown  plotted  in  Figure  (3.10).  Both  curves  show 
considerable  fluctuations  below  5kHz,  due  primarily  to  the  poor  resolution 
of  the  spectrum  analyser  at  low  frequencies.  However,  the  underlying 
shape  of  the  amplitude  ratio  curve  is  similar  to  that  obtained  from  the 
loudspeaker  source  measurements. 

Measurements  of  the  direct  and  transmitted  pulses  were  made 
for  receiver  depths  of  3,  7  and  10  cm  and  the  results  are  shown 
in  Figure  (3.11).  The  puls*  at  a  depth  of  10  cm  shows  considerable 
attenuation  and  is  only  just  discemable  above  the  background 
noise  level. 


Amplitude  Ratio  dB  Phase  Sound  Pressure  Level  dB 
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Figure  3.  10  Anplltude  Ratio  of  Surface  and  Sub-surface 
Spectra. 
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c)  Surface  Mlcrophon 
Figure  3.11  Surface 
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4  GROUND  AS  A  RIGID  POROUS  MEDIUM 

4.1  Introduction 

The  theories  of  spherical  wave  transmission  and  reflection  at  porous 
boundaries  that  have  been  advanced  in  the  previous  chapters  rely  upon 
the  assumption  that  the  ground  can  be  regarded,  effectively,  either 
as  a  locally  reacting  fluid  or  as  a  homogeneous  isotropic  fluid 
sustaining  a  single  (compress ional }  wave  type.  If  the  ground  is 
locally-reacting  then  for  computation  of  the  reflected  field  it  may 
be  characterised  acoustically  by  its  normal  (complex)  surface 
impedance,  and  for  computation  of  the  transmitted  field  the  required 
additional  acoustical  parameter  is  the  propagation  constant  within  the 
effective  fluid.  For  a  model  of  external  reaction  then  computation 
of  the  reflected  field  requires  knowledge  of  both  surface  impedance 
and  propagation  constant. 

A  model  of  the  ground  as  supporting  a  single  compressional  wave  type 
is  consistent  with  propagation  in  a  rigid  porous  medium  in  which  the 
acoustic  disturbance  is  transmitted  only  through  the  pore  fluid  and  is 
attenuated  by  viscous  and  thermal  effects.  This  model  has  been  found 
adequate  to  represent  the  behaviour  of  porous  boundaries  in  room  and 
duct  acoustics  [4.1,  4.2]  and  In  studies  of  outdoor  sound  propagation 
[4.3]. 
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Indirect  estimation  of  the  surface  Impedance  of  any  ground,  that  can 
be  modelled  as  locally  reacting.  Is  possible  from  a  propagation 
measurement  with  a  fixed  source-receiver  geometry  [4. 3, 4. 4].  The 
Impedance  values  result  from  fitting  the  data  on  excess  attenuation 
between  source  and  receiver  after  allowing  for  spherical  divergence 
air  absorption.  Since  the  surface  Impedance  is  complex,  the  fitting 
procedure  requires  simultaneous  variation  of  real  and  Imaginary  parts. 
The  questions  arise  (a)  of  a  suitable  initial  guess  and  (b)  of 
suitable  ranges  and  the  type  of  relationship  between  real  and 
imaginary  parts.  Empirical  formulae  for  Impedance  and  propagation 
constant  In  terms  of  a  single  parameter,  the  specific  flow  resistivity 
of  the  porous  medium  [4.5],  originally  developed  from  acoustical 
measurements  on  glass  fiber  materials, fix  the  frequency  dependence  of 
the  real  and  imaginary  parts  of  the  complex  variables  and  the 
relationship  between  them.  They  have  been  used  to  derive  the 
characteristics  of  ground  surfaces  needed  for  predictions  of 
propagation  near  to  outdoor  ground  surfaces  [4.4,  4.6,  4.7]  In  terms 
of  a  best-fit  flow  resistivity.  Although  successful  in  some  Instances, 
the  empirical  formulae  have  been  shown  to  have  shortcomings  in  this 
application.  These  are  evident  particularly  where  measured  values  of 
flow  resistivity  are  available.  In  these  circumstances  the  best-fit 
value  of  flow  resistivity  has  been  found  to  be  approximately  half  of  the 
mean  measured  value  [4.4],  The  large  scatter  In  the  deduced 
Impedance  values  with  frequency,  also  obscures  whether  or  not  the 
empirical  formulae  predict  the  correct  frequency  dependence.  The 


transmitted  field,  which  is  of  particular  interest  in  this  report, 
depends  upon  the  propagation  constant.  Although  the  empirical 
formulae  include  a  set  for  real  and  imaginary  parts  of  propagation 
constant,  these  have  not  been  tested  for  soils  and  sands  to  the  same 
extent  as  those  for  surface  impedance.  For  these  reasons  the 
remainder  of  this  chapter  is  devoted  to  the  development  of 
theoretically-based  alternatives  to  the  empirical  formulae  for  the 
acoustical  characteristics  of  rigid  air-filled  granular  media  and  to 
a  consequent  examination  of  the  limits  of  applicability  of  the 
empirical  formulae.  Comparison  of  the  theoretical  predictions  and 
measured  data  are  made  for  glass  fibre  materials,  sands  and  soils. 

4.2  The  Modified  -  Raleigh  Model 

4. 2.1. Introduction 

A  comprehensive  review  of  the  available  theories  for  propagation  in 
porous  media  is  given  elsewhere  [4.8].  In  granular  media  the  micro¬ 
structure  is  best  described  in  termsof  grain  shape  size  and  packing 
density,  all  of  which  may  be  determined  accurately  by  non-acoustic 
means.  However  theoretical  approaches  which  specify  the  form  of  the 
rigid  solid  constituent  while  attempting  to  describe  the  acoustic 
field  in  the  fluid-filled  interstices  have  been  constrained  in  each  case 
to  particular  grain  shapes  and  find  particular  difficulties  in 
describing  the  field  near  points  of  grain  contact.  There  is  a  further 
fundamental  problem  with  such  approaches,  of  course,  where  the  rigid 
porous  medium  is  consolidated  such  that  individual  solid  elements 
cannot  be  identified.  A  conceptual  model  first  suggested  by  Rayleigh, 
of  a  rigid  solid  matrix  through  which  run  parallel  Identical  circular 


cylindrical  pores  normal  to  the  surface,  has  been  shown  capable  of 
generalisation  to  many  types  of  porous  media  [4.8].  The  generalisation 
is  achieved  at  the  cost  of  Introducing  an  adjustable  parameter,  the 
oscillatory-flow  pore  shape  factor,  which,  strictly,  is  determinable 
only  from  best-fit  to  acoustical  measurements.  On  the  other  hand, 
a  suggested  semi -empirical  relationship  between  this  parameter  and 
another,  the  steady-flow  pore  shape  factor,  which  is  calculable  for 
certain  known  pore  geometries,  means  that,  in  princiDle,  the  acoustic 
properties  of  a  rigid  fluid-saturated  granular  medium  can  be  deduced 
from  four  non-acoustic  measurements.  These  are  of  flow  resistivity, 
volume  porosity,  grain  shape  and  mean  grain  diameter. 

4.3.  Propagation  in  a  single  pore 

In  this  development,  viscous  and  thermal  effects  are  treated  separately. 
Zwikker  and  Kosten  [4.2]  have  shown  that,  at  least  in  the 
limiting  cases  of  low  and  high  frequencies,  such  Independent 
treatments  give  the  correct  result  for  the  field  within  a 
cylindrical  tube  with  rigid  but  heat-conducting  walls  containing  a 
viscons,  heat  conducting  fluid.  Viscous  effects  may  be  imcorporated 
into  a  complex  (l.e.  frequency-dependent)  density  of  fluid  in  the 
pore  (pc  (u)  )  subjected  to  an  oscillatory  disturbance  with  time- 
dependence  exp(-lut),  (ii  being  the  angular  frequency.  The  result  of 
analysis  may  be  written  [4.2]  as 

Pc  («)  -  P0  [1  -  2  Uc  /f)'1  T  (k^)]*1  (4.1) 

where  T(*  A)  -  J.  («  /TJ  /  J  U  A)  (4.2) 

v  iv  0  C 

JQ( )  and  J,()  being  cylindrical  Bessel  functions  of  zeroth  and  first 
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order,  respectively,  and 

<c  •  a  (w/v)i  (4.3) 

where  a  is  the  pore  radius  and  m  is  the  Kinematic  viscosity  of  the 
fluid  and  pQ  is  its  equilibrium  density. 


In  the  context  of  developing  a  general  theory  for  acoustic  propagation  in 
fluid-saturated  granular  poro  elastic  media,  Biot  [4.10]  considered 
propagation  in  a  paralled-sided  infinite  slit  of  semi-width  b. 

The  expression  for  the  complex  fluid  density  in  such  a  slit 
corresponding  to  equation  (4.1)  may  be  deduced  to  be 

psU)  =  ofl{1  -  [ks  H)*]'1  tanh  [^(-i)*])'1  (4.4) 

where  *  b(w/v)^  (4.5) 

Analysis  of  the  net  flow  of  heat  per  unit  mass  within  a  circular 
cylindrical  tube  containing  a  non  viscous  conducting  fluid,  supposing 

(i)  that  heat  is  transported  in  transverse  (radial)  direction  only  and 

(ii)  that  pressure  is  uniform  over  the  pore  cross  section,  produces 
a  relationship  for  the  complex  (frequency-dependent)  compressibility 
of  the  fluid  within  the  pore,  as  follows  [4,2]; 

Cc(w)  =  (vp^1  [1  ♦  2(r  -  1)T(K*r  *c/i)]  (4.6) 

where  N  is  the  Prandtl  number  and 
pr 


C(w)  =  1  dp 

p  and  p  being  the  time-dependent  density  and  pressure 
similar  analysis  for  net  flow  within  a  parallel-sides 
thickness  2b  produces  [4.11] 


(4.7) 

respectively.  A 
infinite  slit  of 


(4.8) 
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Ce(w)  *  (yPJ'1  (1  + 

3  W  T 

Njr  ks(-D 


Y-1 


tanh  [N*rks(-1)J]) 


The  functions  for  complex  density  for  the  fluid  in  the  circular 
cylindrical  tube  and  the  parallel-sided  slit  have  similar  frequency 
dependence.  The  frequency  dependence  of  the  functions  for  complex 
compressibility  are  similar  also.  The  real  part  of  p(w)  and  the 
imaginary  part  of  C(w)  show  very  little  frequency  dependence  for  a 
wide  range  of  ic.  However  the  variations  of  imaginary  part  ofp  (u,) 
and  the  real  part  of  C(u)  with  k  show  that 


(i)Imps(u)  for  a  slit  of  semi-width  b  a  Imp _ (u> )  for  a  circular 
cylindrical  tube  of  radius  b/n 

and  (ii)  ReC  (w)  for  a  slit  of  semi-width  0  *  ReC  (w)  for  a  circular 

5  C 

cylindrical  tube  of  radius  b/n. 


The  value  of  n  for  which  these  approximations  hold  is  roughly  constant 
and  equal  to  0.5  for  1.8s  *c  s  10.66  (see  Table  4.1).  This  suggests 
that  for  a  pore  of  a  general  shape  but  with  a  constant  characteristic 
dimension  i  along  its  length,  it  is  possible  to  choose  an  equivalent 
radius  a  «  t/n  and  to  evaluate  p  (k)  and  C(k  )  where  k  *  U/n)(u>/v)*  (4.9). 
n  represents  an  oscillatory-flow  pore-shape  factor. 


4.2.3  Extension  to  bulk  medium 

If  the  pores  in  the  bulk  medium  may  be  assumed  identical  and  to  run 
parallel  to  the  surface  then  extension  to  the  bulk  medium  from  the 
expression  for  a  single  pore  Is  straight  forward  and  Introduces  the  volume 
porosity  and  flow  resistivity  in  addition  to  the  single  pore  parameters 


i  and  n.  In  a  real  medium,  for  various  reasons,  including  changing  pore 
orientation  and  cross  section,  the  particle  path  lines  will  deviate  from 
the  surface  normal.  The  changing  pore  orientation  is  referred  to  as  the 
sinuosity  or  tortuosity,  q,  of  the  pore  [4.2,  4.12,  4.13].  The  tortuosity 
factor  q  has  the  particular  value  1/cos  e  for  a  medium  containing  parallel- 
cylindrical  pores  inclined  at  angle  8  to  the  surface  normal.  It  is 
possible  to  incorporate  a  further  factor  (>1)  when  the  pores  have 
cross-section  which  vary  along  their  lengths.  However  this  possibility 
is  not  pursued  further  in  the  present  analysis. 

On  the  assumption  that  all  of  the  pores  are  identical,  the  fluid  velocity 
averaged  over  a  single  pore's  cross  section,  <  v  >,  can  be  related  to 
the  velocity  u,  averaged  over  unit  cross-section  of  the  porous  medium  by 
<  v  >  =  q  u  /a 
where  n  is  the  volume  porosity. 

Consideration  of  the  pressure  gradient  along  a  streamline  in  a 

single  pore  during  laminar  flow  gives  [4.13] 

-  d £  =  q2  4p  u  (4>n) 

dx  0 

♦p,  the  flow  resistance  per  unit  length  in  a  single  pore,  is  given  by 
♦p  =  8us/iZ  (4.12) 

where  u  is  the  dunamic  viscosity  of  the  fluid,  and  s  is  a  steady-flow 
pore-shape  factor,  s  may  be  calculated  from  a  standard  fluid-dynamical 
analysis  for  a  pore  of  any  regular  cross-section.  For  a  circular- 
cylindrical  pore  s  has  the  value  1  and  for  a  parallel  sided  slit 
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s  =  1-5.  A  more  extensive  list  of  values  of  s  corresponding  to  pores 
of  various  regular  cross  sections  is  given  in  ref  [4.13]. 

From  equations  (4.11)  and  (4.12)  the  flow  resistance  per  unit  length,  or 
flow  resistivity  of  the  porous  medium,  which  Is  a  measurable  quantity,  can 
be  calculated  from 

<)>  =  8uq2s/t2Q  (4.13) 

By  eliminating  the  characteristic  dimension  t  between  equations  (4.9) 
and  (4.13)  it  is  possible  to  deduce  that 

k  =  (s*/n)  (8P0q2  *♦)*  (4.14) 

for  pores  of  arbitrary  cross-section. 

The  one-dimensional  equations  of  continuity  and  motion  in  the  rigid 
porous  medium  which  are  sufficiently  general  if  the  medium  is  Isotropic 
may  be  written  [4.2]  as 


and  -  9fi_  =  (q2/tl)  p  (w)  3.u  (4.16) 

a  x  ST 

Combining  these  equations  with  use  of  (4.7)  enables  production 
of  the  compressional  wave  equation  inside  the  rigid  porous  medium. 
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Finally  from  equations  (4.1),  (4.6)  and  (4.17)  the  propagation 
constant  k  within  the  medium  is  given  by 

k2  *  q2(u./co)2  [1  -  2  (k  /T)'1  T  (k/i)  f1 

x  [1  ♦  2  (y  -  1)(M*  k/T)'1  T(N*  k/T)]  (4.18) 

pr  P> 

and  the  relative  surface  impedance  Z/oc  =  1/6  which  is  the 
characteristic  impedance  of  a  semi-infinite  medium  is  given  by 
Z  -  up  (u)/k  (4.19). 

It  is  of  particular  interest  to  note  that  for  low  frequencies  and 
small  values  of  the  characteristic  dimension,  i.e.  small  pores  or  high 
flow  resistivity,  equations  (4.18)  and  (4.19)  reduce  to  relatively 
simple  expressions  by  virtue  of  the  small  argument  approximations  of 
the  Bessel  Functions.  Hence,  at  low  frequencies  and  for  high  flow 
resistivities 

Z/Pac„  *  ks  (1  +  i)  (n2ft6Y(4/(!2Po  sc2)  )i  (4.20) 

4.2.4  Calculations  of  tortuosity 

Tortuosity  has  an  influence  not  only  upon  fluid-flow  in  porous  media 
but  also  upon  their  electrical  conductivity  when  the  pore  space  is 
filled  with  an  electrically  conducting  liquid  and  upon  diffusion  of  water 
vapour  through  the  pore  space.  It  has  been  related  also  to  the  index 
of  refraction  of  fourth  sound  when  the  porous  medium  is  fixed  with 
liquid  helium  [4.14].  In  the  context  of  difusion  of  gases  through 
porous  granular  materials,  the  following  relationship,  known  as 


Bruggeman's  relationship,  has  been  found  to  hold  [4.15], 
q2=ffn'  (4.21) 
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where  n'  is  a  grain  shape  factor  taking  the  value  0.5  for 
spherical  grains.  This  and  other  values  of  n'  are  listed  in  Table 
4.2.  Hence  knowledge  of  porosity  and  grain  shape  should  be  sufficient, 
in  principle,  to  calculate  the  appropriate  tortuosity  value.  It  should 
be  noted  that  for  porosities  near  unity,  the  tortuosity  will  also  be 
near  unity  irrespective  of  particle  shape. 


j 
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4.2.5  Calculations  of  shape  factors 

For  packings  of  spheres,  of  porosities  In  the  range  0.34  <  n  <  0.45 

2 

CARMAN  [4.13]  gives  a  range  of  values  of  2q  s, 

4.5  <  2q2s  <  5.1 

Using  the  Bruggeman  relationship  (4.21)  for  spheres,  the  stated  range 

of  porosities  implies  a  range  of  tortuosities  given  by 

1.49  <  q2  <  1,715. 

2 

Hence  from  the  range  for  2q  s,  it  Is  possible  to  deduce  s  *  1.5  which 
is  the  value  for  parallel-sided  slit  pores. 

In  principle  it  should  be  possible  to  deduce  values  of  s  for  any 
regular  microstructure  from  an  analysis  of  steady  flow  through  the 
known  pore  cross-section,  n,  the  oscillatory-flow  shape  factor  would 
require  computation  from  a  separate  analysis  of  the  microvelocity  field 
for  each  pore  cross-section.  Nevertheless  the  values  of  n  and  s  for 
the  supposed  extremes  of  circular  cylindrical  tube  and  parallel-sided 
slit  suggest  the  possibility  of  the  relationship 

n  -  2  -  s  (4.22) 


I 
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Consequently  for  a  packing  of  spheres,  with  s=1.5,  it  is  possible 
to  expect  that  n  =  0.5. 

For  non-spherical  and  non-uniform  particles  appropriate  values  of  n 
and  s  are  more  difficult  to  specify.  Indirectly  if  the  flow 
resistivity,  porosity  and  mean  grain  shape  are  known  then  it  may  be 
possible  to  deduce  a  value  for  s,  and  hence  n  through  equation  (4.22), 
by  utilising  the  concept  of  hydraulic  radius. 

Hydraulic  radius,  m  is  given  by 

m  =  n/S  (4.23) 

Where  S  is  the  exposed  particle  surface  area  per  unit  volume  of 
material.  A  more  usual  concept  than  flow  resistivity  in  describing 
flow  through  porous  media  is  permeability.  As  long  as  the  pores 
have  fairly  smooth  perimeters,  the  permeability  B  is  related  to  m  by 
B  =  am2  /  (2q2s)  (4.24) 

Flow  resistivity  and  permability  are  related  through 

B  =  n2/*  (4.25) 

Since  m  is  calculated  from  geometrical  considerations  then  knowledge  of 
ifi,  n  and  n'  should  enable  calculation  of  s. 
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4.  3  Comparison  of  rigid  porous  model  with  reflection  measurements 

Although  there  are  many  measurements  available  in  the  literature 

on  the  acoustical  characteristics  of  materials  of  interest  in  room 

acoustics  and  duct  acoustics,  there  is  a  paucity  of  data  for  the  acoustical 

characteristics  of  soils  and  sands.  An  exhaustive  set  of  data,  on 

various  types  of  sand  and  artificial  porous  granular  media 

consisting  of  spherical  lead  shot,  were  obtained  by  Ferrero  and 

Sacerdote  [4.15].  Their  measurement  technique  involved  measurement 

of  the  pressure  standing  wave  (plane  wave  fronts)  formed  in  a 

cylindrical  tube  with  the  sample  forming  one  termination  and  the 

loudspeaker  source  of  pure  tones  forming  the  other.  The  propagation 

constant  and  characteristic  Impedance  were  calculated  from 

measurements  of  surface  Impedance  for  two  different  thickness  of 

sample}  one  thickness  being  twice  the  other.  This  method  was  first 

advocated  by  Pyett  [4.16].  Ferrero  and  Sacerdote  also  measured  the 

porosities  and  flow  resistivities  of  their  samples.  Thus  for  the 

spherical  lead  shot  and  sand  samples  (n’*0.5)  of  all  the  parameters 

are  known,  with  the  exception  of  the  pore-shape  factors  for 

application  of  the  theory  developed  in  earlier  sections  of  this 

2 

chapter.  The  appropriate  value  of  q  for  all  samples  Is  approximately 
1.6,  since  all  a  *  0.4.  Best  agreement  is  obtained  with  the 
£oi  shape  factor  ratio  0.5  rather  than  0.75  or  1.  This  value 

corresponds  approximately  to  n  ■  0.5,  s  »  1.5  l.e.  the  values  for 
parallel-sided  slit  pores  as  expected  for  packings  of  spheres  (see 
section  4.6).  The  measured  data  and  theoretical  predictions  are 
shown  in  Figures  4.1  to  4.5. 
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4.8  Comparison  with  transmission  measurements 

4.8.1  Soils  and  sands 

The  difference  in  sound  pressure  levels  between  a  microphone  at  the 
surface  and  a  probe  microphone  (figure  4.!*)  at  various  depths  in  several 
soils  and  sands  has  be**«  measured  [4.17].  Assumin'?  that  th«  ground  medium 

is  locally  reactino  then  this  difference  should  be  oiven  by 
attenuation  in  dB  =  20  Im(k)  D/loge10 
Im(k)  refers  to  the  imaginary  part  of  k,  i.e.  the  attenuation  constant, 
and  D  is  the  depth  of  the  probe.  Figures  4.6  to  4.9  show  the 
comparison  between  the  measured  data  and  the  predictions  obtained 
through  (i)  equation  (4.19)  (ii)  measured  values  of  flow  resistivity, 
using  an  apparatus  described  elsewhere  [4.18]*  and  (iii)  measured 
values  of  porosity,  deduced  by  weighing  samples  of  known  volume. 

4.8.2  Fibreglass 

Using  the  apparatus  described  in  Chapter  3,  measurements  have  been 
made  of  the  penetration  of  sound  from  an  air-jet  source  into  fibreglass 
50  cm.  thick.  The  probe  receiver  Microphone  was  at  1,4  and  7  cm. 
beneath  the  surface.  The  comparison  between  measured  one-third  octave 
averaged  data  and  predictions  based  upon  measured  flow  resistivity 
(using  apparatus  conforming  to  ASTM)  and  the  assumptions  that 
n  *  q  =  n//T  a  1. 

are  shown  in  Figures  4.10  and  4.11.  Note  that  the  results  on  pulse 
penetration  (Figure  3.17)  are  consistent  with  the  air-jet  results  and 
with  prediction.  In  the  case  of  the  loudspeaker  generated  narrow-band  pulses 
predictions  of  both  the  real  and  Imaginary  part  of  the  transfer  function 
between  surface  and  buried  receivers  (see  sections  2,3  and  chapter  3)  are  in 
good  agreement  with  measurement  (Figures  4.12  to  4.14). 


PROBE  DEPTH/CM 


gp  uoiiDnu»MV 


Frequency  Hz 

Measured  attenuation  in  washed  sand,  o  *  210,000  MKS  Units,  *  0.4  (dry) 
Predictions  use  these  values,  together  with  n’  »  0.5,  n//s“ «  0.75 


Predictions  use  these  values,  together  with  n*  =0.5,  n/ZsT- =  0.75 


tasured  attenuation  in  "institutional^  soil  a  *  ©04,000  MKS  units.  Predictions  use  this  value 
iqether  with  JJ  *  0.4,  n'  «  1.0  and  levelled  values  of  n/fs 


Measured  attenuation  in  Loess,  a  =  650,000  MKS  units.  Predictions  use  th: 
n'  =  1.0  and  labelled  values  of  rt/Ss. 
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5.  GROUND  AS  A  PORO "ELASTIC  MEDIUM  :  RESPONSE  TO  PLANE  WAVES 
5.1  Biot  theory  and  modification  thereto 

For  a  more  complete  description  of  ground  response  to  an  external 
acoustic  disturbance  it  is  necessary  to  model  the  ground  as  a 
poro-elastic  medium. 

we  are  concerned  with  the  acoustic  wave  motion  in  ground,  consisting 
of  solid  particles  and  air  in  the  interconnected  voids.  The  solid 
particles  are  assumed  to  be  in  firm  contact  with  each  other  and 
thus  constitute  a  homogeneous  elastic  frame.  The  theory  of  wave 
motion  in  such  composite  media,  as  derived  by  Biot  [5.1] ,  is  widely 
applicable,  and  predicts  two  types  of  cospressional  waves  and  one 
shear  wave. 

In  this  section,  we  are  interested  in  computing  propagation  constants 
for  the  above  three  types  of  body  waves,  suitable  for  numerical 
calculations  and  for  parametric  studies.  For  this  purpose,  it 
is  convenient  to  express  the  equations  of  motion  in  terms  of  two 
scalar  potentials  and  $2  and  two  vector  potentials  jpi  and  . 

Let  u  and  U  describe  the  displacement  of  the  skeletal  frame  and 
fluid  respectively.  Then  displacement  of  the  fluid  relative  to 
the  frame  is 

w  -  0(u  -  U) 

Then  we  can  express  the  two  vectors  as 
u  «  7fi  +  curl 


w  »  V+2  ♦  curl 
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The  coupled  equations  of  motion  for  tt\e  scalar  potentials  $1  and  $2 
are 

HV2*!  -  CV2*2  *  P*i  -  Pf*2 
CV2$i  -  MV2$2  =  pf$i  -  m$2  -  ^2 

and  another  pair  of  coupled  equations  of  motion  for  the  vector 
potentials  ij^  and 

Gbvzifcj  =  Pii.  -  pfia. 

^  hi  “  Pf  ~  “tz. 


where  H,C  and  M  are  elastic  constants  expressed  in  terms  of  elastic 
moduli  of  the  grain-solid,  fluid  and  frame. 

H  “  <*.  -  V'  J  (D  •  V  +  *b  +  “V3 

C  ”  Kr  (Kr  “  V  '  (D  -  Kb) 

M  »  K*  /  (D  -  K^) 

and 

D  «  K  II  ♦  n<Kr/Kf  -  1)] 

K  is  the  bulk  modulus  of  the  individual  grains,  K  is  the  bulk 
r  1 

modulus  of  the  pore  fluid  ,0^  is  the  shear  modulus  and  1^  is  the 
bulk  modulus  of  the  assemblage  of  particles.  In  Eq.  (5.1),  p 
is  the  average  mass  density,  p^  is  the  fluid  density,  h  is  the 
dynamic  viscosity  coefficient  and  k  is  the  nermeabilitv  of  the  porous 
frame  with  dimensions  {X.]*  . 


(5.1a) 


(5.1b) 


*1 


(5.2) 
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The  parameter  a,  related  to  the  fluid  action,  i*  sometimes  called 
effective  fluid  density  and  is  given  by 

■  “ 

where  ft  is  porosity  and  q  is  tortuosity.  The  parameter  q  (>  1)  accounts 
for  the  fact  that  the  fluid  flow  is  not  in  the  direction  of  macroscopic 
pressure  gradient.  Experimentally  measured  values  of  q  for  soils  and 
sand  have  been  reported  in  literature  [5.2].  However  in  the  present 
study,  we  use  q  as  a  model  parameter  through  the  relation 


where  n* is  a  grain  shape  factor  obtained  from  studies  of  gas  and  solute 
diffusion. 


i 

I 

1 


1 

! 
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In  the  high  frequency  range,  where  the  viaosity  effects  begin  to  be 
of  important,  we  replace  the  ratio  ^  by  n  FU)/k  where  F(X)  is  a 
complex  viscosity  correction  factor,  given  by 

fm  .  l  St  t(a/T) 

r(X)  (  -  1/4)  u  _  2.o  T(xA)/(x/i)T 


and 


X  •  (/a/n)  /8q*  ftkw/u 


(5.3) 
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where  y  is  kinematic  viscosity  (n//s)  is  identified  as  a  single  model 
parameter,  i.e.  shape  factor  ratio.  Note  that  ,  from  equation  (4.1), 
the  complex  density  of  fluid  in  the  pores  of  a  rigid  fluid-saturated  porous 
medium  may  be  written 

p(X)  =  -4FU)  /  (X/iMX/i)] 

We  seek  solutions  of  the  form 

♦  ^  “  A^  exp  (i  (£x  -  uit)] 

and 


$  =  A  exp [I  (4x  -  ut)] 

2  2 

and  Eqs.  (5.1)  provide  a  pair  of  simultaneous  equations  for  Ai  and  A2. 
For  a  non-zero  solution,  we  require  that  the  determinant  should  vanish 
and  hence  we  obtain 


pu>2  -  Hf1 
Pfo>*  -  CV 


Cf-  pfw * 


=  0 


Hi1  -  mu’  +  iuF(X)nA 


The  above  equation  is  quadratic  in  V  and  has  two  distinct  roots 
•  l-B  t  /B1  -  4 AC]  /  2A 


1 


d 


as  long  as  B1  -  4AC  ^  0 
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A  «  C3  -  HM 

1 

B  =  u2{mH  -  2Cpf  +  p  H)  -  itiF  (X)nHA 

C  >=  ul,{pf1  -  mp}  +  o>2 {  ipmF(X)nA)  (5.6)  . 

The  two  roots  of  V  as  given  by  (5.6)  correspond  to  "fast"  and 
"slow"  types  of  compress ional  waves.  The  corresponding  propagation 
constant  for  shear  wave  motion  can  be  obtained  by  assuming 

1 

=  B  exp[i(lx  -  oit)],  «  Bjexplidx  -  uit)]  (5.7)  j 


Then  it  follows  from  the  equations  of  motion  (5.1b) 


pui2  -  G  P  p  w* 


pf<aJ  mci>2  +  iuF(X)Ti/k 


«  0 


II 

) 

(5.8)  -1 


The  phase  velocity  ,w 'Re  l  and  attenuation,  Xm  Z  of  three  types  of 
waves  can  be  easily  computed  from  (5.6)  and  (5.8).  A  FORTRAN  program 
has  been  written  to  compute  propagation  constants  for  two  dilatational 
waves  and  one  shear  wave.  The  complex  frequency  correction  factor 
F ( X )  has  been  computed  using  Bessel  functions  of  coop  lex  arguments 
where  the  calculations  are  performed  in  single  precision. 


II 

II 
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The  pore  parameter  X  as  given  in  Eq.(5.3)  enables  us  to  consider  the 
influence  of  various  parameter  on  phase  velocity  and  attenuation  of 
fast  and  slow  compressional  wave  modes.  For  the  present,  we  restrict 
the  sensitivity  analysis  only  to  the  parameters  relating  to  the  micro 
structure  and  assume  base  values  for  the  eisstic  constants  of  the 
frame  which  are  typical  of  sand. 

A  few  remarks  concerning  the  microstructural  parameters  are  in  order. 
We  have  already  note  that  the  tortuosity  is  dependent  on  porosity 
and  this  gives  us  two  independent  parameters  n‘  and  R  .  Although 
permeability  k  and  shape  factor  ratio  are  related,  in  the  present 
study  we  treat  them  as  another  pair  of  independent  parameters.  Thus 
we  choose  the  four  parameters  Si,  n',  k  andn/Zs  to  describe  the 
microstructure  of  the  ground.  (For  a  more  detailed  discussions  on 
this  see  references  (5.3]). 


Predictions  o*  wave  speeds  and  dispersion  in  air  and  water  saturated  sand 


The  program  described  in  previous  section  has  been  utilized  to  compute 

phase  velocity  and  attenuation  in  water  filled  saturated  sediments  [5.3] 

where  the  values  assumed  for  the  various  parameters  used  in  the 

computation  are  listed  in  Table  1.  Here  we  present,  in  Figures  5.1  to 

5.4  a  sensitivity  analysis  of  dispersion  and  frequency  -dependent  viscous 

attenuation  (fast  wave)  to  porosity  (0.3),  grain  shape  factor  (n*  =  1.25) 
-7  2 

permeability  (10  cm  )  and  shape  factor  ratio  (0.6).  The  values  given 
in  the  brackets  refer  to  base  values. 

The  sensitivity  of  attenuation  to  permeability,  shape  factor  ratio 
and  grain  shape  factor  are  similar  to  those  obtained  for  variations 
in  permeability  and  pore  sisa  parameter  and  for  structure  factor 
respectively  by  Stoll  and  Bryan  [5.4J.  VIbcous  attenuation  appears  to 
be  more  sensitive  to  porosity  as  it  affects  tortuosity  (structure  factor) . 

TABLE  1 


o 


At  tenuat  i  on  C  Neper  s/rrO 


in  grain 


permeability. 


Figure  5.4  Sensitivity  gf  fast  wave  characteristics  in  water-saturated  sand  to  variation  in  shape  factor  ratio 
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In  Fig  5.5(a),  the  prediction  for  the  frequency-dependent  phase  velocity 
and  in  Fig  5.5(b)  the  effect  of  variation  in  frame  elasticity  for  the 
fast-wave  in  air-filled  sand  are  shown.  The  predicted  attenuation  is  less 
than  Q.ldB/m  in  all  cases.  In  Fig  5.6,  we  present  real  and 

imaginary  parts  of  propagation  constants  for  slow  wave  for  the  sake 
of  comparison  with  those  obtained  using  rigid-porous  theory.  This 
shows  that  the  slow  wave  corresponds  to  that  predicted  in  the  modified 
fluid  model  of  a  rigid  porous  medium. 


The  range  of  values  of  oscillatory-flow  pore-shape  factor  (see  section 
4.3)  used  in  these  calculations  (equation  (5.3))  differs  from  that 
quoted  in  Chapter  4  since  it  is  derived  from  comparison  the  forms  of 
F(l)  for  oscillatory  flow  in  a  cylindrical  pore  and  parallel-sided 
slit  rather  than  those  of  p (u) .  The  relevant  parameter  may  be  labelled 
n (viscous)  where 

0.75  4  n (viscous)  4  1 

The  bulk  modulus  of  air,  Kf ,  must  be  complex  to  account  for  thermal 
effects  in  air-filled  porous  media.  The  required  complex  value 
may  be  obtained  from  equation  (4.6)  where  Kf  *=  i/C(w)  and  the  required 
value  for  oscillatory-flow  pore  shape-factor,  n (thermal)  is  such 
that 

0.5  4  n  (thermal)  4:  1 

Relationships  equivalent  to  (4.22)  may  be  invoked  between  the  various 
shape  factors  vis. 

2n (viscous)  *  n (thermal)  +  1 
2n (viscous)  ■  3  -  s 

and 

1  4  s  $  1.5 


where  s  is  the  steady  flow  shape  factor  defined  in  section  4.4. 

For  the  computations  in  Figures  5  *nd  5.6,  «•  have  used  the  base  values 
n(thermal)  «  0.5  and  n (viscous)  »  0.75.  The  remaining  parameter 

values  <u;e  given  in  Table  II. 

Note  tnat  a  single  value  of  oscillatory-flow  pore-shape  factor 
n (viscous)  may  be  used  if  C(»)  is  calculated  from 

CM  -  (YPJ_>  12  -  Y  ♦  2(y  -  1)  W^X/i)  /TW^X/i)) 

where  y  is  the  ratio  of  specific  heats,  is  the  Prandtl  Number  and 


po  is  atmospheric  pressure. 


TABLE  II 


Parameters  Required  By  Biot  Theory 
(As  Modified  By  Attenborough/Rayleigh) 


Parameter 

Grain  Bulk  Modulus 

Grain  Mass  Density 

Frame  Bulk  Modulus 

Frame  Rigidity  Modulus 

Fluid  Bulk  Modulus  (isothermal) 

Dynamic  Viscosity  of  Fluid 

Porosity 

Permeability 

Grain  shape  factor 

Shape  factor  ratio  (n//a) 


“Sand"  Value 


3.6  x  !0U  dyn/cm2 


2.6B 


g  cm 


8.27  x  ^0^  dyn/cm2 
8.27  x  10^  dyn/cm2 


.00184 
0.4 

3.76  x  10“6  cm2 
0.5 
0.6 


dyn-sec/cm 


Rmm  Velocity  (a/aoc) 


1150.00 


1000.00  4 - - - ' — ' — *  "*'*  ‘  ' - ' - ' — ' — 1  ’  '  '  ' 

0.01  froqMWKy  -  kHz  l.C 

Figure  5.5(a)  Predicted  fast-wave  phase  velocity  for  air-filled 
sand  with  characteristics  given  by  Table  II. 


I  - - - - - - - — . . * 

0.01  fTMMWDT  -  kHz 


Figure  5.5(b)  Variation  of  fast-wave  speed  characteristic  in 
air-filled  sands  with  labelled  values  of  fraae 
elastic  constants  I 


1.00 


0.01  frapwncy  -  kHz  i.Of 

5.6  Predicted  real  and  imaginary  parts  of  slow- wave  propagation 
constant  in  air-filled  sand  with  characteristics  given  in 
Table  II. 


As  has  been  shown,  the  values  predicted  for  the  fast  wave  propagation 
constant  In  air-saturated  sand  are  very  sensitive  to  the  value  taken  by 
the  frame  rigidly  (shear)  modulus  (assumed  equal  to  the  bulk  modulus) . 

The  shear  modulus  for  soils  is  affected  by  various  soil  parameters  and  the 
state  of  stress.  The  important  parameters  include  [5.7]  effective 
octahedral  normal  stress,  void  ratio,  grain  characteristics,  shape,  size, 
grading,  mineralogy,  amplitude  strain,  ambient  stress  history,  secondary 
time-dependent  effects,  degree  of  water  saturation,  octahedral  shear  stress, 
frequency  of  vibration,  soil  structure  and  temperature.  This  complexity 
implies  difficulty  with  selecting  values  typical  of  any  given  soil  type. 
Consequently  three  values  have  been  chosen  which  indicate  a  typical  range 
and  enable  a  study  of  the  sensitivity  of  propagation  characteristics  to 
shear  modulus  in  air-filled  soils.  The  values  chosen  [5.6]  are  set  out  below :- 

_2 

Soil  type  Shear  Modulus  (G)  dyn  cm 

O 

Santa  Barbara  clay  1.38  x  10 

AGSCO  No. 2  5.17  x  10® 

9 

Ottawa  sand  1.86  x  10 

Figures  5.7  to  5.10  show  the  relative  speeds  and  attenuations  of  the  three 

possible  wave  types  predicted  by  assuming  equal  shear  and  bulk  moduli  and 

using  the  base  values  of  other  parameters  given  in  Table  II  with  the 

exceptions  r. '  -  1  and  shape  factor  ratio  »  0.75.  Fioure  5.7 

shows  that,  for  a  relatively  high  value  of  rigidity  modulus,  the  'slow' 

wave  is  clearly  the  slowest  up  to  1kHz,  the  fast  wave  having  a  frequency 

independent  value  somewhat  greater  than  the  velocity  of  sound  in  air 

and  similar  to  the  value  predicted  for  dry  sand.  From  Figures  5.7  and 

5.8  it  is  clear  that  a  small  decrease  in  rigidity  modulus  causes  the  shear 

wave  to  become  faster  than  the  slow  wave  near  1kHz.  Note  that  both  of 

the  rigidity  modulus  values  assumed  in  generating  Figures  5.7  and  5.8  produce 

fast  wave  speeds  somewhat  less  than  tha  speed  of  sound  in  sir. 


Further  decrease  in  the  rigidity  modulus  to  the  value  assumed  in  Figure  5.10 
produces  the  interesting  situation  where  the  'slow'  and  'fast*  waves  are 
predicted  to  exchange  roles  above  100  Hz.  It  should  be  noted  that  still 
lower  values  of  rigidity  modulus  are  possible  for  soils  and  these  will 
produce  'slow'  and  'fast'  wave  exchange  at  still  lower  frequencies.  Finally 
addition  of  a  complex  component  of  the  shear  modulus  is  shown, in  Figure  5.11 
to  produce,  as  expected,  much  higher  attenuations  for  the  'fast'  and  shear 


waves . 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUACAU  OF  STANOAWS  -  ”  A 


Figure  5.7  Predicted  propagation  constants  for  dry  soil  of  rigidity  modulus 
..  -  fast  wave  _ slow  wave  _  _  .  shear  wave 


Attenuation  <dB) 


Figure  S.8  Predicted  propagation  constants  for  dry  soil  of  rigidity  modulus 


Figure  5.10  Predicted  propagation  constants  for  dry  soil  of  rigidity  nodulus 
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5.3  Predictions  of  surface  solid  particle  velocity  resulting  from  plane 
wave  incidence 

The  ratio  of  solid  particle  velocity  to  the  pressure  on  the  surface 
has  been  used  to  describe  the  ground  response.  To  facilitate 
computations  of  this,  we  make  use  of  a  different  but  equivalent 
formulation  to  that  presented  in  section  1.  The  present  formulation 
enables  us  to  express  both  solid  and  fluid  displacement  in  terms  of 
two  scalar  potentials  4^  and  and  two  vector  potentials  4  and  4 
We  write 

u  -  grad  4j  ♦  curl  4i 

U  «  grad  4^  ♦  curl  **  (5.9) 

and  they  obey  the  following  coupled  differential  equations 

7‘(*i  ♦",»  *  a£(W  +bra)  -V  1 

V2(Q4,  ♦  *4  )  -  —t  (p ,4  )  -  »l»)  1^(4  -4  )  (5.10) 

\  I  1  2 

and 

MVi4i.|^(pi4i)tbr(l)|l<^-4i, 

°  ’  S'  <P«!l)  •  br(X)  Ji  \  "V  «*•«> 

where  b  •  n/fc,  and  p^  are  densities  of  solid  and  fluid  respectively. 

The  elastic  perasmters  p,  q,  it  can  be  related  to  (5.5). 

bulk  mod<.' lus  of  grains  ,  K{  (bulk  modulus  of  fluid),  1^  (bulk  modulus 

of  aggregate) 


J 
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N  is  the  shear  modulus  of  the  solid  and  composite 


<1  -  fl)  [II  ■  HI  -  ^  ]  l  ♦  H  «  1^ 


+  4, 


Q  « 


i  -  a  -  «■  n  xr 

F" 


H  -  n  -  jn  k 

tT 


l" 


U  -  fl  -  ^ 


+  K 


Kr  Kf 


a2  x 


[l  -  a  -  +  n  xrl 


Prom  the  computational  point  of  view,  the  present  formulation  does 
not  pose  any  additional  problem  as  the  previous  program  can  be 
still  used  and  Eqs.(5.l2)  provide  for  the  translation  between 
the  two  notations. 


(5.12) 


J 

I 

1 

1 


3! 

31 

f 

fl 
fl 


For  comparison  with  maasuramsnts  it  is  convenient  to  calculate 
the  square  of  the  ratio  of  surface  particle  velocity  to  pressure 
henceforth  called  the  intensity  ratio.  These  calculations  are 
described  for  obliquely- incident  plane  waves  on  a  semi-infinite 
poro-elastic  medium  and  normal-incidence  on  a  rigidly-backed 
poro-elastic  layer  in  Appendix  4. 

Figures  5,)2  and  S.13  show  plots  of  normal -incidence  intensitv  ratio 
against  frequency  for  the  parameter  values  given  in  Table  II. 

These  are  expected  to  be  representative  of  sand.  It  should 
be  noted  that  the  number  of  layer  resonances  increases  with 
increasing  depth. 

The  variation  of  intensity  ratio  with  angle  of  incidence,  on  a 
semi-infinite  layer  of  dry  sand,  measured  away  from  the  normal 
and  with  frequency  is  shown  in  Figure  5.14.  The  ratio  decreases 
with  increasing  angle  of  incidence  and,  for  any  given  angle,  the 
ratio  decreases  with  increasing  frequency. 


Intensity  ratio  X  1-0E-08 
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THE  SURFACE  MOTION  OF  A  PORO- ELASTIC  HALF-SPACE  DUE  TO  AN  INCIDENT 
SPHERICAL  SOUND  WAVE. 


Introduction 


while  the  treatment  of  the  ground  surface  as  a  rigid  porous  medium 
gives  an  adequate  description  of  the  reflected  sound  wave  and  the 
transmitted  slow  wave  (consisting  of  air  motion  within  the  pores) , 
any  attempt  to  determine  the  seismic  response  must  take  into  account 
the  finite  rigidity  and  density  of  the  solid  medium. 

Here  we  take  account  of  frame  elasticity  by  using  the  theory  of 
propagation  in  poroelastic  media  outlined  in  the  previous  chapter. 

This  has  been  developed  extensively  by  Biot  in  a  number  of  publications 
In  these  he  does  not  always  use  the  same  scheme  of  notation,  and  we  opt 
for  the  system  given  by  him  in  1961  [6.1].  This  gives  the  most 
compact  expressions  in  subsequent  derivations,  and  many  of  the  material 
parameters  are  defined  in  such  a  way  as  to  be  identical  to  those  of  the 
corresponding  rigid  porous  material. 

The  problem  of  spherical  wave  incidence  on  a  poro-elastic  half  space 
is  approached  in  the  same  manner  as  that  for  the  more  well-known 
cases  of  fluid-fluid,  fluid-solid  and  solid-solid  interfaces,  which 
which  are  described  by,  for  example,  Ewing  et  al  [6.2].  The 
axisynmetric  displacement  potentials  are  expressed  in  terms  of 


inverse  Hankel  transforms,  the  transformed  potentials  being  determined  by 
the  boundary  conditions  at  the  interface.  In  this  case  there  are  four  potentials, 
corresponding  to  the  reflected  congressional  wave  in  the  upper  medium,  the 
transmitted  slow  and  fast  conpressional  waves,  and  the  transmitted  sheer  wave. 

The  coaiplexity  of  the  resulting  integrals  prohibits  a  thorough  analysis  using 
the  saddle  point  method  as  in  the  fluid-fluid  case  used  to  model  a  rigid  porous 
ground  surface.  However,  a  simpler,  less  rigorous  approach  enables  us  to 
obtain  asymptotic  expressions  for  the  various  wave  types  at  grazing  incidence, 
and  it  is  also  straight  forward  to  identify  the  propagation  paths  and  evaluate 
times  of  arrival  for  transient  sources. 

Some  simplification  is  possible  by  considering  the  limiting  case  when  the 
bulk  modulus  and  density  of  the  saturating  fluid  are  much  smaller  than  the 
elastic  modulus  and  density  of  the  aolid.  The  propagation  constants  of  the 
fast  and  shear  wave  are  unaffected  by  the  presence  of  the  fluid,  while  that 
of  the  slow  wave  la  the  same  as  for  the  rigid  porous  case.  This  'light 
fluid'  limit  should  be  applicable  for  air-filled  ground  surfaces. 

Peng  and  Johnson  (6.3,  6.4]  have  recently  considered  the  problem  of  a 
fluid/porous  solid  interface.  Their  treatment,  otherwise  comprehensive, 
gives  results  for  the  high  frequency  limit  when  the  propagation  constants 
for  the  various  wave  types  are  all  real  and  non-dispersive .  This  allows 
them  to  define  'true'  surface  waves  with  a  velocity  less  than  any  of  the 
body  waves,  and  to  i^lemant  the  Cagniard  method  whan  evaluating  the 
response  to  a  pulsed  source.  However,  the  simplest  possible  treatment 
of  ground  surfaces  must  admit  a  complex  alow  wavs  propagation  constant,  and 
as  a  result  ws  era  forcsd  to  give  s  mors  general  analysis  of  surface  wave 


contributions. 


( 
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6.2.  Displacement  potential  in  the  upper  medium 

i 

The  vector  of  fluid  displacement  in  the  upper  medium  is  defined  in  terms 
of  the  potential  t|>: 

v  “  Vip  (6.1) 

The  acoustic  pressure  is  given  by 

p  -  -KaZ.-»  ■  -KaV2*  <6-2> 

where  Ka  is  the  bulk  modulus 
ip  obeys  the  wave  equation 

V2*  -  <6-3> 

everywhere  in  the  upper  medium  except  at  the  source. 
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6.3  Displacement  potentials  in  the  lower  medium 


If  the  displacement  vector  of  the  solid  matrix  is  u,  then  following 
Biot  [6.1]  we  represent  the  displacement  of  fluid  in  the  pores  relative 
to  the  solid  by  the  vector 


w  »  0(0  -  u) 


where  0  is  the  porosity  and  U  is  the  actual  average  fluid  displacement. 
Hence  v  is  the  relative  volume  displacement  of  the  fluid  over  unit 
area  of  the  bulk  medium. 

Because  the  solid  can  support  shear  stress  the  displacement  are  not 
in  general  irxotatlonal.  They  can  be  detained  in  terms  of  scalar  and 
vector  potentials  aa  follows: 

w  •  ♦  V  x 

u  •  ♦  2.  *  x. 

The  coupled  equations  of  propagation  for  dilatatlonal  waves  are 

■  ft  “>♦, ♦  w 
♦  "V  ■  ft  (p,*,  ♦  +  fw- 

and  for  shear  waves 


(6.4) 


(6.5) 

(6.6) 

5.7a) 

(6.7b) 


and  a 


is  a  parameter  which  relates  changes  in  pore  volume  to  overall 
changes  In  volume. 
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6.4  Boundary  condition#  at  the  interface. 

The  definitions  of  fluid  and  solid  displacement  used  here  lead  to 
particularly  simple  expressions  for  the  boundary  conditions. 

Assuming  open  pores  at  the  interface,  the  boundary  condition  for 
continuity  of  fluid  displacement  is 

n.  (u  +  w)  »  n. 

where  n  is  a  unit  vector  in  the  vertical  direction.  In  terme  of  the 
displacement  potentials  this  becomes 

n.V  ♦  O  ♦  n.J.  «  (Xj  ♦  ■  n .V*  (6.9) 

Biot's  notion  of  'effective  stress'  is  particularly  appropriate  for 
the  application  of  stress  boundary  conditions.  The  effective  stress, 
t^,  is  that  part  of  the  total  stress  within  the  porous  msdium  (acting 
on  both  fluid  and  solid)  in  excess  of  the  local  fluid  pressure.  The 
components  of  effective  stress  acting  in  the  plane  of  the  interface 
must  therefore  be  equal  to  sero.  Follvoing  Biot's  definitions  and 
using  Cartesian  coordinates  x,y,s  with  the  s  axis  in  the  vertical 
normal  direction,  we  have 

Ss  "  2w*s  *  "  2«  “  ■*>•  ♦  Hit  -  o)  C 


0 


(6.10) 
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t* 

yz 


0 


(6.111 


t* 

xz 


0 


(6.12) 


where  e 

z 


3  2 


e  =  V^u 

5  =  -Z.JL 


Finally,  the  pressure  in  the  pores  must  be  equated  to  that  in  the  upper 
medium,  giving 

p  ♦  one  -  m;  =  0  (6.13) 


Me  later  express  the  stress  boundary  conditions  in  terms  of  the  potentials 
when  transforming  to  axieyometric  coordinates. 
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The  complex  density  P  is  identical  to  that  for  the  equivalent  rigid 
2 

porous  material.  The  reduction  to -this  case  can  be  achieved  simply  by 
setting  the  solid  displacement  potential  4^  equal  to  zero  in  equation 
(6.15b).  This  gives 


Prom  equation  (4.16)  we  can  deduce  the  complex  density  for  a  rigid  porous 
material  and  we  have 


J 

I 

.1 


Axisvmnetric  boundary  conditions 


Since  the  bracketted  tern  must  approach  zero  as  r  ■+• 


this  boundary 


condition  means  that  it  must  be  zero  everywhere,  that  is 


+  V2* 


0 


After  substituting  the  appropriate  potentials,  equation  (6.13)  becomes 

mV2*  ♦  aMV2*  -  1C  V2*  *  0 
2  s  a 

He  shall  ultimately  use  the  four  boundary  conditions  (6.21)  -  (6.24) 
to  determine  the  four  potentials  4,  4^ ,  4^,  4^. 


(6.23) 


(6.24) 


I 
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6.8  General  solutions 

Assuming  a  point  source  of  sound  at  a  height  z0  above  the  interface,  the 
solution  to  (6.14?  may  be  written 


ikjR 

4-  *vr*z> 


where  R’  =  r'  ♦  (z  -  l,)1  ,  and  4^  represents  the  wave  reflected  from  the 
interface. 

The  general  solution  in  axisymmetric  coordinates  may  be  written  as  an  inverse 
Ranke 1  transform: 


♦  ^(r.zJ  ■ 


1  (K)e"VlZ  J  (Kr)XdX 

l  * 


(6.25) 


where  v  *  (K1  -  k')1,  n  v  >  0,  t  >0,  Im(k  )  >  0.  and  4  is  the 
ill  1  l 

transformed  potential  for  *  -  0. 


The  total  field  above  the  interface,  using  the  integral  representation 
for  the  point  source,  is 


B-vi | »-*. | 

vl 


♦  ?!• 


-vxs 


J.  (Kr)KdK 


(6.26) 


B 
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k'j  is  the  propagation  constant  for  the  equivalent  rigid  porous  material, 
and  k',  is  the  propagation  constant  for  congressional  waves  in  the  bulk 
medium  assuming  that  relative  motion  between  fluid  and  solid  is  constrained 
in  some  way. 


Equation  (6.32)  Is  a  quadratic  which  may  be  solved  in  the  usual  way 
to  give  two  possible  values  of  a. 


If  the  second  term  in  (6.32)  is  negligible  then  we  must  have 


(6.35) 


For  our  application  we  do  in  fact  expect  that  (6.32)  will  give  two 
wavenumbers  corresponding  to  distinct  wave  types,  and  also  that  kj  <  kj 

Accordingly,  we  will  call  one  solution  k, ,  and  this  is  approximately 
equal  to  kj  and  associated  primarily  with  motion  of  air  in  the  pores 
known  as  the  'slow  wavs'.  The  second  solution  we  call  k, ,  and  this  is 
the  wavenumber  of  the  compressions 1  'fast  wave'  determined  chiefly  by 
the  properties  of  the  solid  part  of  the  medium. 


I 
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In  equations  (6.29)  and  (6.30)  the  two  potentials  must  be  in  a  fixed 
proportion  for  each  allowed  wavenumber.  When  a  »  k,  we  use  (6.29) 
to  obtain 


tlTpj  -  k?aM 
H(k*  -  kj*) 


and  when  a  =  k,  substitution  into  (6.30)  gives 


iL  =  -  m2P1 

i,  mk;1  -  kj) 


Equations  (6.15k),  (g.l5b)  therefore  have  two  independent  solutions. 


For  the  slow  wave 


(k)eV2  i  J0(Kr)KdK  ) 


and  for  the  fast  wave 


4,  =  ♦,'(*)  eV,Z  Jc(Kr)KdK  ) 


(6.39) 


*2  -  *».♦! 


v2  -  (K*  -  k*)3 

v,  -  <**  -  ^ 
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V  i*  that  part  of  tha  tranafonaad  potential  $2  associated  purely 
with  the  slow  wave.  Similarly  $'  Is  the  part  of  ^  associated 
with  the  fast  wave. 


In  general  we  may  write 


(*)ev,z  ♦  a, ,  ♦J,eV*  *)  J.(Kr)KdX 


(i,,eV*  *  +  a,1i,,eVj*)  J,  (Kr)Kdk 


Finally,  the  general  axisyenetrlc  solution  for  the  shear  wave 
motion  described  by  equations  (6.19)  is 
«• 

9  -  f.$(lt)eV**  J, (kr)KdK 


where  z<0,  v  »  (Ka  -  kfc) 
* 


(6.40) 


(6.41) 


(6.42) 


v 


J 


r 
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Simultaneous  equations  for  tho  transformed  potentials 

The  general  solution  given  by  equations  (6.26),  (6.40)  -  (6.42) 
for  2  =  0  «u:e  now  substituted  into  the  boundary  conditions 
(6.21)  -  (6.24).  This  leads  to  four  linear  equations  in  the 
transformed  potentials  which,  when  satisfied,  also  satisfy  the 
boundary  conditions. 


i 


Equation  (6.21)  is  satisfied  if 


v  i  +  v  (1  +  «!»)$,’  +  u  (1  +  a,  ,)$,*  -  K1  (1  - 
1*1  I  J  P2  a 


Equation  (6.22)  is  satisfied  if 

(2ua,,v^  -  kj  [a,,  (H  -  2jj  -  oM)  -  M(1  -a)]} 
+  {2jiuJ  -  kj  £  H  -  2p  -  aM  -  a,,M(l  -  a)]} 

-  2yv  K*  |  =0 

Equation  (6.23)  is  satisfied  if 

2«, .^i,1  ♦  2v,$  -  (2K>  -  k;)4  -  <J 


Equation  (6.24)  is  satisfied  if 


+  kjM(l  +  aa,  ,)£  +  kjM(a  + 
e-VlZ0 

V i 


(6.46) 


where  we  have  used  k1  «  uj2  Pi 


with  some  rearrangement  we  now  write  equations  (6.43)  -  (6.46) 
in  matrix  form: 

-VlZO 

(A]  {*}  «  {S} 

This  is  shown  in  full  in  equation  (6.47) 

Using  Cramer's  rule  the  transformed  potentials  are 


1  A . 


Ai  e 


-vizo 


A»  V! 


«!? 

h 


A a. 
Ao 


Aa. 

Ao 


♦  .a  A» 

’*  r 


~ 

Vi 

B~VjZo 


(6.48) 


(6.49) 


(6.50) 


(6.51) 


where  Ao  is  the  determinant  of  the  matrix  (A),  and  Aj,  Aj ,  etc  are  the 
determinants  formed  by  replacing  the  appropriate  column  with  the 
S  . 


vector 
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Solutions  in  integral  fora 

Substitution  of  the  transformed  potentials  into  the  general  solutions 
(6.25),  (6.40)  -  (6.42)  gives 


Ai  e 
Ao 
o 


-vj (z  +  ro) 


Vl 


-Jo  (Kr)KdK 


.  -visj  +  VjZ 

_i  5 - Jo  (Kr)  XdK 

Ao  Vi 


♦  a. 


A.  e 


-VlZo  +  VjZ 


Ao  Vi 

-VjZo  ♦  VjZ 


-Jo (Kr)KdK 


• _ 

j  A,  V! 


-Jo (Kr)KdX 


♦  «. 


Ai  a'v»z°  *  v»z 

Ao 


Jo (Kr)KdK 


♦  •  - 

* 


At,  e 


-VjZO  ♦  v*z 


Ao  vi 


Jo (Kr)KdK 


(6.52) 


(6.53) 


(6.54) 


(6.55) 


where  z  >  0  for  and  z  <  0  for  the  other  potentials. 


These  Integrals  are  typical  of  those  found  in  other  two-media  probl, 
and  we  may  consider  them  qualitatively  in  terms  of  saddle  -  point 
integrals,  branch  cut  integrals  and  poles  in  the  complex  k  plane. 
The  saddle  points  occur  at  the  stationary  points  of  the  oscillating 


integrand,  and  for  large  Kr 


of  these  correspond  to  the  ray  paths 


governed  by  the  laws  of  reflection  and  refraction.  The  remainder 
consist  of  Inhomogeneous  waves  (see  Brekhovskikh  [6.5])  for  which 
there  is  no  definite  path,  but  which  can  be  interpreted  in  terms  of 
Snell's  law  by  considering  complex  angles  of  incidence  and  refraction. 

The  branch  cut  integrals  correspond  to  rays  which  are  in  part  tied  to 
the  interface.  Some  of  the  remaining  parts  of  these  paths  are  indefinite 
since  they  consist  of  inhomogeneous  waves,  but  again  interpretation 
is  possible  by  allowing  complex,  as  well  as  real,  critical  angles. 

when  a  pole  of  the  integrand  lies  between  the  contour  of  integration  and 
the  real  axis  there  will  be  a  surface  wave  contribution.  Nearby  poles  may 
also  affect  the  saddle-point  integrals  as  in  the  rigid  porouB  case. 

The  eventual  contour  of  integration  may  not  always  consist  of  all 
possible  saddle-point  and  branch  cut  integrals.  If  the  source  and/or 
receiver  are  sufficiently  distant  from  the  interface,  a  single  saddle 
path  may  enclose  more  than  one  branch  point,  thereby  excluding  any 
separate  contributions  from  the  associated  branch  cut  integrals. 
Physically,  this  appears  to  eliminate  any  path  which  consists  partly  of 
inhomogeneous  waves  and  whose  time  of  travel  is  less  than  the 
corresponding  wholly  homogeneous  path.  At  grazing  angles,  however, 
each  saddle  path  encloses  a  single  branch  point,  and  all  possible 
arrivals  have  to  be  considered. 
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These  considerations  are  suassarised  diagramatically  in  Fig.  (6.1) 
in  a  manner  closely  following  that  given  by  Ewing  et  al  [6. 2)  for 
a  solid-solid  interface.  For  this  example  we  assume  k^  >  k^  >  k^  >  k^. 
The  wavy  lines  denote  inhomogeneous  waves,  and  in  general  two-part  paths 
from  branch  cut  integrals.  The  potential  +1  is  associated  purely 
with  the  slow  wave  and  the  potential  with  the  fast  wave,  so  that 
in  equations  (6.531,(6.54)  we  could  write 

*7  =  +  •».♦»»  ♦>  •  +  «>>♦> 

In  addition  to  these  arrivals,  of  course,  there  may  be  surface  wave 
contributions. 

It  can  be  seen  from  the  diagram  that  when  the  receiver  is  on  the 
interface,  the  four  paths  for  each  branch  point  associated  with  the 
four  different  wave  potentials  coalesce  into  a  single  path.  This 
reduces  the  total  number  of  paths,  excluding  surface  waves,  from 
sixteen  to  four,  and  this  is  a  useful  simplification  when  we  consider 
motion  of  the  ground  surface.  Expressions  for  this  are  derived  in 
Appendix  5. 
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CHAPTER  7 

CONCLUDING  REMARKS 

7.1  Although  we  have  been  able  to  predict  qualitative  features  of 
some  existing  data  on  acoustic-to-seismic  coupling,  in  particular, 
that  used  to  compute  intensity  ratios  (Chapter  5) ,  the  requirement 
for  design  of  unattended  seismic  and  acoustic  sensor  logic  remains 
that  good  quantitative  agreement  with  these  ratios  should  be  obtained 
and  that  there  should  be  theoretical  agreement  with  other  data  on 
acoustic-to-seismic  coupling  17.1).  These  data  include  the  results 
of  tests  using  impulsive  acoustic  sources  above  and  below  the  surface 
and  sensors  at  depths  of  up  to  1.5  m. 

7.2  The  main  obstacles  to  such  requirements  at  this  stage  are  lack  of 
information  on  the  frame  elastic  constants  required  to  characterise 
poro-elastic  soils  (Chapter  5)  and  the  complex- layered  and  possibly 
anistropic  structure  of  real  grounds. 

7 . 3  The  analysis  has  shown  that  surface  particle  velocity  prediction 

is  particularly  sensitive  to  the  value  chosen  for  the  frame  rigidity 
modulus. 

7.4  The  model  developed  in  this  work  has  been  brought  to  the  stage  where 
numerical  predictions  of  coupling  from  above-ground  point  acoustic 
sources  to  sub-surface  geophones  are  possible  in  certain  well-defined 
situations  vie.  near-surface  sources  and  air-saturated,  isotropic, 
porous  ground  of  infinite  depth.  Alter  nr  .-ively ,  if  the  acoustic 
source  is  sufficiently  far  removed  from  the  ground  surface  for  the 
incident  sound  waves  to  be  treated  as  plane,  then  not  only  can 
predictions  be  made  of  the  coupling  to  the  surface  of  semi-infinite 
ground  but  also  to  ground  in  the  form  of  a  hud-backed  porous  layer. 


7.5  Extension  of  this  model  to  allow  for  anistropy  and  for  more 
complicated  ground  structures  and  continuous  sources  will  be 
straight  forward  albeit  sufficiently  algebraically  cumbersome  to 
warrant  a  semi-numerical  approach,  such  as  that  adopted  by  the 
Fast  Field  Program  and  its  successors  in  Underwater  Acoustics. 

7.6  Further  extension  of  the  model  to  allow  for  impulsive  rather 

than  continuous  sources,  should  also  be  possible  using  a  numerical 
approach  as  suggested  recently  by  Feng  and  Johnson  [7.2] . 

7.7  The  model  may  also  be  extended  to  include  the  effects  of  partial 
water-saturation  of  poro-elastic  soils  along  the  lines  suggested 
by  Brutsaert  [7.31. 

7.8  Finally,  it  should  be  noted  that  the  theories  developed  in  this 
report  are  applicable  to  seismic  disturbance  of  poro-elastic  ground 

by  direct  mechanical  contact.  Consequently  it  is  possible  to  envisage 
prediction  and  analysis  of  seismic  refraction  results  in  a  way  that 
differs  from  established  methods.  Furthermore  the  possibility  arises 
of  the  prediction  of  ground  vibration  due  to  road  and  rail  vehicles 
in  terms  of  basic  fundamental  attributes  of  intervening  terrain. 


J 
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WEYL  -  VAN  DER  POL  FORMULA  --  CONFIRMATION 

The  method  of  steepest  descents  has  been  employed  extensively  to  obtain 
a  computationally  feasible  solution  for  the  field  from  a  point  source 
of  sound  in  air  above  a  locally  reacting  ground  plane.  Various 
mathematical  approaches  which  arrive  at  the  above  solution  under  various 
approximations  are  listed  in  reference  [1]. 

The  important  case  of  near  grazing  incidence  requires  additional  care  while 
handling  the  steepest  descent  integrals  and  this  is  chiefly  achieved  by  a 
method  known  in  the  literature  12]  as  subtraction  of  the  pole.  This  method 
often  leads  to  laborious  calculations.  A  simple  alternative  solution  can 
be  obtained  by  an  appropriate  application  Of  the  method  of  stationery  phase.  In 
general, the  method  of  stationary  phase  yields  useful  approximation  as  long 
as  the  integrand  is  a  slowly  varying  function  of  integration  variables  in  the 
neighbourhood  of  the  stationary  point.  After  a  suitable  choice  of  integration 
variables,  the  integrals  are  approximated  in  the  usual  way  and  the  final  results 
for  the  reflected  field  confirm  those  obtained  previously. 

The  physical  problem  under  consideration  can  be  formulated  in  terms  of  the  following 
boundary  value  problem.  For  convenience,  we  choose  a  cartesian  coordinate  system 
The  acoustic  velocity  potential  satisfies  Helmholtz  wave  equation  in  both  the 
media  (z  <  0  and  <  0)  which  are  characterised  by  their  densities  Pi  and  P2 
and  propagation  constants  kj  and  <z ,  respectively.  Thus  we  have. 


distances  by  the  method  of  stationary  phase  and  we  obtain 


! 


R(9)  e~ilciR* 
4  it  r2 


where  R2 
2 


»  x2  ♦  y2  +  lz  t  z  )2 


R(e> 


cos  e  -  m 

COS  0  +  M 


J~2  2 

/n  -  sin  e 
/n^  -  sin2  0 


and  M  ■  p^/pj  and  n  *»  .  Hie  above  asymptotic  approximation 

is  of  limited  use  and  breaks  down  when  0  is  near  s/2  and  this  is 
due  to  the  fact  that  the  integrand  in  (A1.3)  no  longer  remains  a 
slowly  varying  function,  but  will  oscillate  rapidly  between  -1  to  +1 
In  this  case,  <0  *  */2J ,  the  contribution  from  the  stationary  point 
has  to  be  considered  more  carefully.  In  order  to  facilitate  this 
computation,  we  rewrite  the  integral  in  (A1.3)  as 


e~iklR2 
4*  R„ 


(1) 


*  * 


where 


(1) 


4w 


lax 


da 


P-X 
1  2 


Wi^xV 


-iX,  (r-z  J+lBy 

A  O 


(A1.4) 


The  integration  is  performed  successively,  first  with  respect  to 
(from  now  on  abreviated  as  w.r.t.)  8  and  then  w.r.t.o.  The  integral 
w.r.t.  6  is  evaluated  after  deforming  the  contour  into  the  cosqplex  6 
plane  and  then  approximating  the  resulting  integral  by  the  method  of 
stationary  phase. 


This  is  accomplished  by  first  observing  the  branch  point  and  pole 
singularities  of  the  integrand  in  the  complex  8  plane.  In  the 
present  case  the  integrand  has  branch  points  at 
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The  deformed  contour  consists  pf  two  circular  arcs  of  large 


radius  going  round  the  branch  cuts.  The  contribution  from  the 


circular  arcs  is  zero.  The  contribution  from  the  branch  point 

,  ■>  h 

~{k  -  a)  can  be  shown  to  be 

2 


0  (exp{-l (k  2  -  a2) **  {***  )-  -  i  tie  2  -  a2)**  y>) 

1  o  2 


and  will  be  neglected.  Thus  we  are  left  with  the  following  integral 


-iX  v  eosfl  +i  Rv  sin8 

«  1  1  .1  1  1 

dS 


1*2 


IX  (p2Yp1X2) 


(A1.5) 


where  we  have  let  y  «  r^  sin  9^,  z  +  zq  •  r^  cos  9^  in  A1.5  and 
C  is  the  contour  round  the.  branch  point  -  (r^2  -  a2)*1.  This 
integral  is  approximately  evaluated  by  the  method  of  stationary 
phase.  The  exponent  in  A 1.5  has  a  stationary  point  at 
B  -  -  Sk2  -  a2  sin  9^.  When  9  v  w/2,  the  situation  we  are 
interested  in,  we  observe  that  the  stationary  point  is  in  close 


proximity  of  the  branch  point 


-  j^c  2  -  a2. 


Now  the  main  contribution  to  the  integral  arises  from  the  neighbourhood 
of  stationary  point,  in  that  part  of  the  integrand  which  is  slowly 
varying.  The  case  where  the  integrand  is  not  slowly  varying  has  to 
be  considered  separately.  In  order  to  do  this,  we  first  rewrite  (A1.5)as 

/—  ■ 


-S  f  R  (a, 6  J  -  1}  {  /k  2  -  a2  cos  0  +  k  a  }  x 

v  i  x  11 
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-i X  (z+z  )+i8y 
e  1  ° _ 

X1 ( *i+*iA) 


d6 


c 


(A1 .6) 


2 

where  A  *=  M  /-? — i-  and 

1-M' 


P  X  P  v 
R  (a  ,9  )  =  2  1  "  12 

v  j  - 

P  V  +  P  y 
2  *1  1  *2 


At  6 


-  /<  2-«2 


sin  8. 


Now  we  write 

X  »  •/<  2  -  a2  -  82  »  (/k  2  -  a2  +  8)  (Xe  2  -  a2  -  8) 

X  l  1  1 

In  the  neighbourhood  of  the  stationary  point,  the  second  factor  is  slowly 
varying  and  we  put  8  •  -  Xr^2  -  a2  (as  8^  *  w/2) . 

On  the  remaining  Integrand,  we  put 

(  £J~7,  .iH  .  Ji  •l'/* ' 

where  x  is  a  real  parameter  and  varies  from  to  +  “  on  C. 

After  some  algebra,  we  have  from  (At. 6) 

{  V“'V  -1  )  (  co*  *!  ♦  -  °2’  A) 


x  exp  {  -  iy^  -  a2  (1  +  •>  cot^a^l  }  x  X  Qr^sin  *j| 


(A1 


where 


I(u)  - 


e~t  u 
t  +b 


dt 


and 


b  «  *j  e  iW^4  f{x  2  -  a2)  cot  9 ,  x  A  (x  2  -  a2)  '} 


hr 


l 


1  l  i 


(Al. 


The  Integral  defined  by  I(u)  can  be  expressed  in  closed  font  involving 
error  function  complement  and  we  can  write  it  as 

I(u)  «  -  /§»  exp  (-iw/4)  ,P(b  exp  (i  s/4)  n^uj 


Thus  we  have  from  (A1.4)  and  (A1.7), 


4ir 


^iux-lyl^2  -  a  2  (1+S  cot 4  9^) 


-S 


«  (coee^x^x^  -a2,  *}  t(  2^  air,  ^ 


Ida  (‘i 


(X) 
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2 

In  the  present  case  of  Interest  y(l  +  *i  cot  0^)  *  y  *  since  cos  0^ 
is  small.  The  final  evaluation  of  ii  ^  is  carried  out  by  a  further 
application  of  method  of  stationery  phase.  Now  the  integral  1(2^  sinSj) 
does  not  contribute  to  the  exponential  and  the  stationary  point  of 

ax  "  x\^  *  ~  “2  occurs  at  a  -  -  Kj*/(x2  +  r2)  .  The  relevant  terms  at 
this  point  are  given  by 


K  t 
1  1 

R 


since  sin  81  *  1 


K  COS  8 
1 


cos 


K 

1 


y 


R  sin  8^ 


b 


0.5  exp  (-  i  w/ 4)  iA  +  cos  9} 


From  (A1.3),  (A1.4)  and  (A1.9),  we  finally  obtain  for  the  reflected  field 
in  z  *  O  as 


-i<  R 
e  1  1 
4ir  R^ 


-iK  R 
«  1  2 
4  if  rr 


1  *  2Ly  (R(01  -1)  F(t)  } 


(A1.10) 


where  y  -  /  *jic  R 

1  2 


(A+  cos  9) 


A  »  M 


/~l 

V  1  -  H2 


*  8  wh«n  )n|2  »  1  and  Im I  «  1 


and  F(a)  -  exp  (ia  ) 


m. 

|  exp  (-it' 


)  at 


After  soee  simple  algebra  the  above  solution  tor  the  reflected  field  is 
easily  shown  to  be  equivalent  to  those  previously  obtained. 
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APPENDIX  2 


ON  THE  ACOUSTIC  SURFACE  WAVE  NEAR  A  RIGID  POROUS  BOUNDARY 

A  spherical  wave  incident  upon  an  interface  of  two  media,  in  particular, 
air  and  porous  ground  gives  rise  to  reflected  and  refracted  waves  subjected 
to  Snell's  law.  Under  certain  circumstances,  the  interface  can  support 
a  wave  which  is  exponentially  damped  in  the  normal  direction  and  whose 
amplitude  decays  as  r  ^  This  wave  is  termed  as  a  surface  wave. 

In  this  section,  we  clarify  certain  mathematical  details  which  are  necessary 
for  the  existence  of  a  surface  wave.  we  also  provide  an  alternative 
interpretation  for  the  acoustic  surface  wave  in  terms  of  discrete  spectrum 
of  a  suitable  partial  differential  operator. 

We  consider  an  incident  spherical  sound  wave 

*1*1 

a. _ 

4» 

with  time  dependence  of  the  form  exp  (-iv>t) ,  incident  upon  a  locally 
reacting  boundary.  For  the  above  case,  surface  waves  are  conveniently 
analysed  by  expressing  the  reflected  field  as  a  Fourier  integral. 

The  pole (s)  of  the  reflection  coefficient  (in  Fourier  domain)  are 
closely  related  to  surface  wave(s) .  The  computation  of  the  reflected 
field  crucially  depends  upon  the  evaluation  of  the  following  integral 
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II 

I 


♦  <H<r>  * 


P 

2 

2u 


-"vlH 


p  V  +  p  V 
2  1  12 


J.(KT)  K  dK 


(A2.1) 


where  v  ■  »KJ  -  k 

l>2 


and  p  are  deneitles  of  the  media  z  >  0  and 

l  2 


z  <  0  respectively  and  H  «  z  +  z„.  The  square  root  functions  are  defined 

in  such  a  way  that  -  k  2  =  i  k ,  , ,  as  K  -e>  0.* 

12  1,2 
* 

Due  to  the  presence  of  brtmch  point  singularities  in  the  integrand  of  (a2 . 1 ) 

it  is  important  to  define  the  branch  cuts  before  analysing  the  existence 

of  poles.  The  choice  of  branch  cuts  is  arbitrary  as  long  as  Re  v  >0 

1.2 

to  ensure  the  convergence  of  the  integral  in(A2,l).  Essentially  two  types 
of  branch  cuts  are  used  in  the  literature  viz,  vertical  and  hyperbolic  [1) 
branch  cuts.  Here  we  consider  the  case  of  vertical  branch  cuts  and  the  analysis 
for  hyperbolic  branch  cuts  follows  on  similar  lines. 


The  vertical  branch  cuts  are  defined  by  He  K  *  t  Re  kj,  Im  K  >  0  and 

Re  K  »  t  k2<  Im  K  >  0.  Corresponding  to  each  square  root  function,  there 

are  two  choices  for  the  sign  of  Re  and  these  are  shown  in  table  *2 . 1 . 

HOLE  A2. 1 


Re  vx 

Re  v 

2 

I 

♦ 

♦ 

II 

- 

+ 

III 

♦  • 

- 

IV 

*k  are  complex  constants  with  non  zero  Imaginary  points 


1.2 
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We  write  v.  «  i  v  I  exp  (i  arg(v  )],  for  i  =  1,2  and  the  sign  of  Re  v 
1  i  i  i 

depends  on  the  variations  in  arg  (v^) .  we  consider  these  variations  by 
considering  a  small  neighbourhood  around  the  branch  points  kj  and  k2 

TABLE  A2 . 2 


e 

Arg  v  (K) 

l 

e 

Arg  (K) 

-3x/2 

-3w/4  ♦  Xl/2 

-3x/2 

-3»/4  ♦  Xj 

•IT  —  X 

A1 

-»/2 

-  Y 
*2 

-x/2 

-»/2  ♦  Xl 

-»  +  X 

-x/2  +  X 

2 

2 

-X1 

0 

-X2 

0 

X1 

X 

1 

X 

2 

?2 

x/2 

x/4  ♦  Xj 

x/2 

x/4  +  xi 

where  Xi 

is  arg  kj . which  is  seiall  as 

kj  fs  essentially  real 

and  x  is  arg  k  . 

Clearly  on  sheet  I, 


and 


-»/2  +  <  arg  (K)  <  x/4  +  x^ 


-x/2  +  x  <  arg  v  (K)  <  x/4  +  X  • 

>  2  2 


Now  we  consider  the  possible  existence  of  poles  on  the  upper  Riemann  sheet. 
It  is  easily  seen  that  the  poles  of  the  integrand  in  (A2.1)  are  given  by  the 
equation 

P  v  ♦  p  v  ■  0 
2  1  1  2 

i.e. 


K  ■  i  ki  I <1  -  Ml  n2)  /  U  -  M*)l* 

where  M  -  p  /  a  and  n  =  k  /k  .  In  order  to  locate  these  poles  in  the 
12  2  1 

first  quadrant,  we  define  the  complex  number 

X„  -  kj  [<1  -  W2n2)  /  U  -  M2)]* 


(A. 2. 2) 
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2 1 

For  most  of  the  grounds  of  practical  interest  | M  j  <<  1  and  we  approximate 

(1  -  K2)*  1.  Bence  K„  -  kj  (1  -  M2n2)J 

We  also  write  the  complex  number  nM  •  | nM |  exp  (ii(i) .  It  is  readily  shown  that 

K0  «  kj  fl  ♦  |nM|4  -  2  !n«|2  cos  2*]*  (A 

and 

arg(K0)  ■  arg(kj)  +  Jarg  <1  -  n2M2) 


-  +  Itan 


Inal  2  sin  2V 
|na|2  cos 2}  -  1 


(A2.4) 


Since  Xj>0  and  K„  >  |kj|,  the  pole  lies  to  the  right  of  the  branch  point 
kj  as  shown  on  Fig  A2.1.  in  order  to  verify  that  the  equation  for  the  pole 
is  really  satisfied  at  K„,  we  compute  Vi(K)  and  vj (K)  at  K„ . 

Vi<K„)  «  t  ix^na  and 
arg{v1(K„)}*  ±  w/2  +  arg  (kj)  +  arg  (nM) 
where  the  choice  of  sign  is  fixed  depending  on  the  Riemann  sheet.  Now 
arg  (kj)  is  a  small  positive  quantity  and  arg  (na)  Si-n/4  for  a  porous 
material  at  low  frequencies  and  for  K0  situated  as  shown  in  Fig.  A2 . 1 . 

Xj  S.  ar9  ^(KJ  s  »/<  +  Xl 

on  the  upper  Riemann  sheet  and  hence  Vi(Ks)  >  +  ik^nM. 

Similarly, 

v2  (K0)  =  ±  ikjn 
2 

where  |n|  is  assumed  large)  end 

arg  V2<KS)  «  ±  ir/2  +  arg  (kj)  +  arg  (n) 

Now,  for  K„  lying  to  the  left  of  the  branch  point  k^,  we  have 

-ir/2  +  x2  <  «rg  V2OC.)  <  0 

Since  arg  (k2 )  >  arg  (kj),  we  have  on  sheet  1  , 

Vj(K.)  *  -  ikjn 

Thus  it  is  evident  that  on  sheet  I,  the  equation  for  the  pole  is  satisfied 
and  hence  K.  is  a  pole  on  sheet  I. 


I 
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Now  the  evaluation  of  the  integral  la  accomplished  by  first  extending  the 

range  of  integration  fro®  0  to  •  to  to  +•  ,  and  deforcing  the  contour  into 

upper  half  plane  of  cut  K  -  plane  which  encloses  the  pole  at  K..  The  branch 

cut  integrals  are  approx  lasted  using  the  saddle  point  nethod  and  the  residue 

term  represents  the  surface  wave  _ 

Aik  exp(-ik  HA)  .  .  /  _ 

residue  at  K,  ■  - = -  8.  (k.r  /  1-A  ) 

(1-M  )  4* 


where 

A  *  Ml  (n*  -  1)  /  (1  -  M*)]* 


For  the  sake  of  completeness  we  quote  the  final  expression  for  the  reflection 
field  in  2  >  o. 


♦ref 


where 


4wRi 


inr  \r  (e*>  * 


2  HI 


*  B  cost).) 


k1R^B  +  cos8.)* 
+  residue  at  the  pole. 


_  A** _ _ 

k,R2 


ocr,  ;  V 


N  =  )IR"  (0.)  +  R'  <0J  COt  eo  ) 

R(8„)  «  (cose.  -  8)  /  (cos  8.  +  0  ) 
6  »  Mn 


At  this  stage  it  is  worth  mentioning  that  the  above  solution  is  of  limited 
practical  use.  For  a  solution  which  is  uniformly  valid  for  all  angles  of 
incidence  a  modified  solution  is  used  as  described  in  previous  section  (see 
Bq.(Al.lO)  ),  or  by  a  siailiar  method  [Al.l).  In  this  case,  the  surface  wave 
term  appears  implicitly  through  error  function  complement,  which  can  be  seen 
more  clearly  for  large  argtaents  of  the  function  (see  Bq.  52  in  (Al.l)). 


J 
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A  simple  interpretation  for  the  surface  wave  can  be  given  in  terms  of  the 
discrete  spectrum  of  a  suitable  partial  differential  operator.  In  the  present 
case  it  suffices  to  consider  a  two  dimensional  model  for  ease  df  illustration, 
where  x  and  z  coordinates  denote  the  horizontal  seperation  and  vertical 
distance  from  the  origin  respectively. 

we  are  required  to  find  a  function  $(x,z)  such  that 

4xJ{  +  Oc'  -  L)4  -«  i(x)  6 (z  -  z 0 )  (A2.5) 

where  L$  =  ,  over  the  half  plane  -<»  <  x  <  «■>,  0  <  z  <  <■>,  such  that 

♦z (x,0)  +  a*  (x,0)  =  0  (A2.6) 

where  a  is  a  paramenter  related  to  Impedance  and  Re  a  >  0.  Also  we 
require  the  field  to  be  outgoing  at  infinity. 


This  operator  espression  can  be  evaluated,  once  we  know  the  spectral 
representation  for  L  under  the  above  boundary  condition  .  This  can  be 

obtained  easily  by  the  well  known  methods  [A. 4]  and  we  have 


6(z  -  z„)  *  2a  exp(-a  (z  +  so))  [1  -  H  (Rea)) 
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Now  from  A2.7  and  A2.8 


-♦  <x,l)  ' 


-1  a  exp  (l^V  +  ?  Ul)  exp[-a(z  +  zj] 
/k1  +  o* 


exp 


(ip(z  -  z.)  +  i*'*1  -  pJ  (x 1 1 


dp 


lip(z  ♦  z.)  ♦  i/k5  -  p5  | x | ]  -£~zJ£  -P—  - 

p  +  la  /k1  -  pJ 

In  Eq.  the  fleet  ten  arises  from  a  pole  of  the  Green's  function  and 

In  the  literature,  1s  referred  to  as  constituting  the  discrete  spectrum  of 
the  operator  L.  the  remaining  two  integrals,  after  some  algebra,  can  be 
easily  seen  to  be  contributing  to  the  continuous  spectrum.  In  fact,  they 
can  be  shown  to  be  direct  (first  integral  in  (A2.9))*nd  reflected  (second 
integral  in  A2.9  waves.  In  other  words  we  notice  that  the  branch  point 
singularities  of  Green's  function  give  rise  to  space  waves,  while  pole 
singularities  give  rise  to  surface  wave (s) . 


(A2.9) 
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APPENDIX  3 

THE  ACOUSTIC  MODELLING  FACILITY 

Introduction 

This  Appendix  introduces  the  .basic  rules  of  acoustic  modelling 
and  discusses  the  essential  requirements  of  the  modelling 
instrumentation.  It  concludes  with  a  description  of  a  series 
of  model  experiments,  conducted  in  a  small  anechoic  chamber, 
to  measure  the  amplitude  ratio  of  the  direct  and  transmitted 
sound  fields  above  a  porous  half-space  and  the  delay  between 
the  arrival  of  the  direct  and  transmitted  sound  waves. 

A3. 1  Basic  modelling  rules 

In  order  to  model  an  acoustic  sound  path,  two  conditions  must 
be  fulfilled.  The  first  is  that  of  geometric  similarity  and 
the  second  kinematic  similarity.  Geometric  similarity  ensures 
that  the  ratio  of  the  wavelength  of  the  sound  source  to  the 
dimensions  of  any  bodies  or  reflecting  surfaces  remains  constant 
and  hence  diffraction  effects  will  be  replicated  in  the  model. 

Kinematic  similarity  is  required  due  to  the  fact  that  air  is 
used  as  the  propagation  medium  in  both  the  full  scale  and  model 
environments  and  hence  the  velocity  of  sound  is  a  constant.  This 
imposes  a  time  scale  factor,  as  can  be  shown  by  considering  the 
dimensionless  quantity  characterising  reference  tines  (the 
Thomson  Number)  as  follows: 


where  T,V  and  L  are  time,  velocity  and  length  respectively  and 
the  subscripts  m  ando  refer  to  the  original  and  model  propagation 
paths.  Since  Vm  =  VQ  then 


Tm 


Lm 

IT 


and  the  time  scale  factor  equals  the  length  scale  factor.  Thus 
maintaining  a  constant  sound  velocity  requires  that  the  model 
frequencies  be  increased  with  respect  to  those  of  the  original 
and  that  all  time  periods  bv  scaled  downward. 


In  general  the  scale  factor  derived  for  the  condition  of  geometric 
similarity  must  also  be  applied  to  the  acoustic  properties  of 
the  materials  used  to  reconstruct  the  surfaces  and  obstructions 
in  the  model ;  so  that  the  impedance  at  the  model  frequencies 
equals  that  at  the  full  scale  frequencies. 


The  prime  objective  of  the  model  experiments  conducted  during 
this  study  was  to  evaluate  theoretical  predictions  of  the 
coupling  between  the  airborne  and  seismic  waves,  rather  than  to 
obtain  quantitative  data  of  a  precise  ground  surface.  In  view 
of  this  it  was  felt  unnecessary  to  employ  'scaled'  models  and  so  no 
attempt  was  made  to  scale  the  acoustic  properties  of  the  model 
ground  surface. 
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A3. 2  Requirements  of  a  model  receiver 

The  basic  requirements  of  a  model  receiver  are  small  size,  high 
sensitivity  and  wideband  frequency  response.  Size  is  important 
for  two  reasons  (a)  because  directivity  is  related  to  the  ratio 
of  diaphragm  diameter  to  wavelength  and  hence  a  small  microphone 
better  approximates  an  omni-directional  receiver  and  (b)  because 
the  reduction  in  diaphragm  mass  improves  the  high  frequency  response 
of  the  microphone.  The  major  disadvantage  of  a  small  diameter  is 
that  the  reduction  in  surface  area  reduces  the  microphone's 
sensitivity. 

The  receivers  used  for  the  model  experiments  were  commercial 
1/4  inch  diameter  condenser  microphones  manufactured  by  Bruel 
and  Kjaer  (B  i  K)  of  Denmark.  The  sensitivity  and  bandwidth 
of  typical  microphones  are  shown  in  Table  A3. 1 


Table  A3. T  Microphone  sensitivity  and  bandwidth 


Microphone 

Nominal 

Sensensitivity 

Bandwidth 

type 

diameter 

mV/Pa 

kHz 

B  &  K  4136 

1/4" 

4 

IQ 

6  (  K  4138 

1/8" 

1 

140 

Figures  A3. 1  -  A3.4  shew  the  directivity  of  these  micmnhones  at 
various  frequencies,  measured  both  with  and  without  a  protection 
grid  [A3. 1 ].  Although  the  1/8"  displays  a  wider  bandwidth  and 
superior  directivity,  its  low  sensitivity  limited  the  maximum 
length  of  the  model  propagation  path  to  less  than  1  metre. 
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In  view  of  this  the  1/4“  microphone  was  adopted  as  the  standard 
receiver  and  the  1/8"  microphone  was  used  primarily  to  test 
the  performance  of  the  1/4"  at  short  separation  distances.  The 
1/8"  microphone  was  also  used  to  measure  the  directivity  of  the 
sound  sources. 

Requirements  of  a  model  sound  source 

The  sound  source  should  be  omni-directional,  have  a  wideband  frequency 
response  and  generate  adequate  signal  power. 

Model  sound  sources  fall  into  two  categories,  continuous  and 
impulsive.  The  impulsive  source,  such  as  an  electrical  spark, 
generates  a  short  energy  pulse  of  high  acoustic  intensity  and 
wide  bandwidth  which  propagates  along  an  infinite  number  of  paths. 
Each  path  can  be  identified  by  the  time  delay  between  the  energy 
discharge  and  the  arrival  of  the  pulse  at  the  receiver,  hence  there 
is  no  need  for  special  anechoic  measurement  environments. 

The  major  limitation  of  the  impulsive  source  is  the  possibility 
of  non-linear  acoustic  propagation  near  the  vicinity  of  the 
source,  which  imposes  a  minimum  length  of  propagation  path. 
Furthermore,  they  require  sophisticated  measurement  and  analysis 
systems  to  cope  with  the  problems  of  direct  and  reflected  pulse 
overlap  [A3. 2]. 

Continuous  sources  can  be  handled  with  relatively  simple  measurement 
systems,  but  since  the  energy  is  radiated  continuously,  it  is  not 
possible  to  identify  individual  propagation  paths.  This  means  that 
all  measurements  must  be  conducted  within  an  anechoic  environment. 


1 
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4  Model  sound  sources 

Two  sources  were  selected  for  the  model  experiments.  The  first 
was  the  air  jet  described  by  Delany  et  al  [A3. 3]  and  shown  in 
Figure  A3. 5.  Measurements  of  the  directivity  in  both  the  vertical 
and  horizontal  planes  indicate  that  this  design  is  a  good 
approximation  to  an  omni-directional  source  at  frequencies  up  to 
100  kHz.  One  particularly  important  feature  of  this  design  is  its 
excellent  long  term  stability,  so  that  over  a  period  of  six 
months  the  one  third  octave  spectrum  varied  by  less  than  0.5  dB. 

The  useful  bandwidth  of  this  source  is  1-100  kHz,  as  dictated 
by  the  signal  to  noise  ratio  of  the  measurement  system  and  the 
low  frequency  performance  of  the  anechoic  chamber.  Figure  A3. 6 
shows  the  free-field  spectrum  measured  at  a  distance  of  100  cm. 

The  second  model  sound  source  was  a  spark  discharge  device  designed 
at  the  Open  University.  (A  complete  description  of  the  design 
and  operating  characteristics  is  given  in  reference  A3.£) 

The  spark  operates  at  very  low  energy  levels  and  so  non-linear 
propagation  effects  are  virtually  eliminated.  A  further 
advantage  of  the  low  energy  is  the  fact  that  electromagnetic 
radiation  is  kept  below  levels  that  can  be  detected  by  the 
condenser  microphones.  A  typical  spark  source  spectrum  is 
shown  in  Figure  A3. 7. 
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A3. 5  Modelling  instrumentation 
A3. 5.1  Dynamic  range 

Acoustic  modelling  is  very  demanding  of  instrumentation  systems. 

The  low  sensitivity  of  small  microphones  requires  that  pre-amplifiers 
and  main  amplifiers  have  very  low  electrical  noise  floors  over  a  very 
wide  frequency  range. 

A  typical  microphone  amplifier  chain  consists  of  a  pre-amplifier, 
to  convert  the  high  output  impedance  of  the  microphone  to  a 
low  impedance  suitable  for  driving  long  cables,  and  a  measuring 
amplifier  to  provide  the  polarization  voltage  and  the  gain. 

An  estimate  of  the  dynamic  range  of  the  instrumentation  system  can 
be  obtained  by  calculating  the  noise  floor  of  the  amplifiers 
and  subtracting  this  value  from  the  Input  signal  level.  The 
following  figures  are  taken  from  the  relevant  manufacturers' 
specifications. 

(a)  Pre-amplifier  noise 

The  inherent  noise  of  the  pre-amplifier  is  determined  by  the 
capacitance  of  the  microphone  cartridge.  For  a  1/4" 
microphone  the  total  noise  over  the  range  20-20  000  Hz 
is  quoted  as  less  than  120  microvolts. 
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(b)  Main  amplifier  noise 

The  noise  floor  of  the  main  amplifier  is  determined  by  the 
thermal  noise  of  the  electronic  components.  The  signal  to 
noise  ratio  for  a  100  mV  input  is  stated  to  be  better  than 
100  dB,  which  implies  a  noise  floor  of  1  microvolt  over  the 
range  20-20  000  Hz. 

Assuming  that  the  main  amplifier  noise  can  be  neglected, 
compared  to  the  pre-amplifier  noise,  then  the  equivalent 
sound  pressure  level  of  the  noise  floor  is  63.5  dB. 

If  non-linear  propagation  effects  are  to  be  avoided,  then 
the  maximum  source  sound  pressure  level  must  be  kept  below 
140  dB.  Subtracting  the  noise  level  from  the  signal  level 
gives  a  dynamic  measurement  range  of  approximately  75  dB. 

A3. 5. 2  Atmospheric  absorption 

The  source-receiver  dynamic  range  is  further  degraded  by  atmospheric 
absorption.  It  is  well  known  that  high  frequency  sound  is 
attenuated  as  it  propagates  through  the  air  due  to  a  combination 
of  thermal,  viscous  and  molecular  relaxation  effects.  Oelany  [A3. 4] 
has  made  an  extensive  study  of  the  attenuation  mechanisms  and 
derived  a  set  of  equations  that  enable  corrections  to  be  made 
based  upon  measurements  of  the  ambient  temperature,  air  pressure 
and  relative  humidity. 
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The  results  of  calculations  for  typical  indoor  values  of  these 
parameters  shows  that  the  attenuation  can  exceed  5  dB/metre 
at  100  kHz. 

A3. 5. 3  Oigital  processing  equipment 

The  majority  of  the  signal  processing  tasks  undertaken  during  this 
study  were  executed  on  a  commercial  digital  spectrum  analyser, 
a  Nicolet  660B.  This  is  a  dual  channel  analyser  performing 
a  Fourier  Transform  on  1024  data  points  to  produce  a  400  line 
spectrum  over  a  pre-determined  measurement  bandwidth.  Alternatively, 
the  instrument  can  produce  a  one  third-octave  band  analysis  over  its 
full  measurement  bandwidth  of  100  kHz. 

The  analyser  has  five  basic  operating  modes,  but  the  only  one 
used  for  this  study,  was  the  dual  channel  mode.  If  the  two  input 
signals  are  represented  by  A  and  B  respectively,  then  the  analyser 
will  calculate  the  cross-spectrum  from  the  expression 

G(f)AB  =  A(f)  .  B(f ) 

and  the  transfer  function  from  the  expression 
H(f)  =  G(f)AB/A(f) 

The  Nicolet  is  also  capable  of  calculating  the  coherence  function, 
which  for  linear  systems  can  be  interpreted  as  the  fractional 
proportion  of  the  mean  square  value  of  the  output  signal  B(f)  which 
is  contributed  by  the  input  signal  A(f).  The  coherence  function,  g  (f) 
can  be  defined  in  terms  of  the  cross-spectrum  as 

i 
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i6<f»»n'* 

I  A(f ) | . |B(f )  | 


Measurement  of  the  coherence  function  proved  a  useful  means  of 
estimating  the  influence  of  near  surface  turbulence  upon  grazing 
incidence  propagation. 


Figure  A3. 5  Air  Jet  Noise  Source. 
(After  Delaney  et  al) 
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APPENDIX  4 


CALCULATION  OF  SURFACE  INTENSITY  RATIOS 


Semi-infinite  medium 


In  this  section  we  are  interested  in  obtaining  an  expression  for  the 
square  of  the  ratio,  solid  particle  velocity  to  pressure  at  the  surface, 
where  the  ground  is  considered  to  be  of  infinite  depth.  We  consider 
a  plane  wave  incident  upon  the  interface  at  an  angle  6, away  from  the 


We  note  that  in  a  poro-elastic  half  space,  there  will  in  general  be 
six  types  of  waves.  These  consists  of  two  dilatational  and  one  shear 
wave  in  both  the  solid  and  fluid  phases.  In  terms  of  scalar  potentials, 
we  write,  for  the  solid  phase 

4'  =  a'  exp  [-iuit  -  k'y  sine'  -  Mx  cos0')l 

is  ^ 


*”  =  a”  exp  [-iuit  -  k*y  sin0"  -  k’x  cos0")  J 
1  s 


i|/  =  aY  exp  [-iuit  -  )fy  sintf  -  l^x  cos$  )) 

i  s 


and  for  the  fluid  phase, 


(A4.2) 


where  k ' ,  k"  and  1^  are  propogation  constant  as  obtained  in  section 
S.l  t',  t"  and  i*  are  the  following  ratios  of  amplitudes 
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and  the  relevant  components  of  solid  strain  are 


,  4  /3u  3u  \ 

3ux,  e  *  3uy ,  e  »  1_  I _ x  -  _x. 

ST  **  3T  **  2  lay  } 


and  finally  the  fluid  dilatation 


3u  3u 

E  =  _ X_  +  _ 

3x  3y 


He 


also  have  the  stress-strain  relations  for  an  isotropic  medium 


o  =  (P  -  anu^  *  2Wtij  + 


i3 


os. j  -  <Q“kk  +  H  V  4ij 


and  inplicit  sunmation  over  repeated  indices  is  assumed. 


The  normal  solid  particle  velocity,  V  is  given  by 


3u 

V  =  x 
s  3t 


and  from  Eq  (5.16)  ,  we  have  for  V  at  the  surface  x  -  0, 


iw  lag'  (  i  k'  cos  8'  ♦  ag"  (i  k"  cos  6") 

+  agr  (  i  kr  sin  8  r) J 
x  exp  (i  k  sin  8  y  -  i  wtj 


Wherein,  we  have  made  uae  of  the  relation. 


k  sin  8  -  k*  sin  8'  -  k"  sin  8-  •  kr  sin  9r 

o  o 


(A4.5) 


(A4.6) 


(A4.7) 


(A4.8) 


(A4.9) 


(A4.10) 
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Also  at  the  surface  the  normal  stress  In  the  frame  must  equal 
stress  In  the  fluid  as  both  are  acted  upon  by  the  outer  air 
presure  p.  Thus 


XX  =  £  »  -p 

(i  -  fl)  n 


From  the  stress-strain  relations  <A4.7)  we  can  express 


-Kf  l|_^e  ♦  e]  j 


and  using  (A4.5)  and(A4.6)  ,  we  get  for  the  pressure  of  the  surface 


P  =  -K. 


r  t  -  a)  (-a  ,  *k'2-a  1  +{-as‘T'k,; 

1  v  a  •  8 


x  exp  [iX  y  sin  e  -  iut) 
o  o 


(A4.ll) 


From  (A4.9)  arid  (A4.ll)  we  obtain  the  ratio  of  solid  particle  velocity 
to  the  pressure  at  the  surface  i.e. 


V  «  U  -  12) li l cob  0’  ♦  I'  k_[_  cos  8'  ♦  I’,  KV  sin  I)  ) 
—  1  k"  . r 

p  _ £ _ 


*rrr 

r/1  -  a  *  t’1 

,♦  r  /t  -  a  * 

n 

a  1 

1  a  J 

'j 

(hi.  I/I 


r  -  a  "  k" 

1  a 


_  r  .  r 

V,1 


The  constants  and  are  as  yet  unknown  and  can  be  expressed  in 
terns  of  known  quant ites  with  the  help  of  the  following  conditions 
at  the  surface  x  *  0.  These  are: 

i)  vanishing  of  tangential  stress  at  x  «  0  i.e.  e^  =  0 

ii)  stress  in  the  frame  equals  stress  in  the  fluid  i.e. 


These  two  conditions  provide  two  equations  for  T  and  T  and  after 

1  2 

some  algebra  we  have  from  (A4.5)  and  (A4.7), 
sin  2  6'Tj  -  cos  2  8r  r2  «  -  sin  2  6' 

and 

(2  nsi  cos2e-  ♦  t*  o  -  (i  -  n>Q]>  Tj  +  mi  sin  2er  r2 
-  2  NO  cos2  0r  -  {Afl  -  (1  -  SI) Q}  (A4.13) 

which  are  solved  to  give 

r!  *  -IMP  tan  2  9r  sin  2  8*  ♦  2Hfl  cos2  9'  ♦  {AO  -  U  -  P)Q  }) 

(Nil  tan  2  0r  sin  2  8"  +  2NP  cos2  0”  +  { All  -  (1  -  P)Q  }  ) 

and 

r2  =  r  sin  2  8"  ♦  sin  2  8 '  (A4.14) 

cos  2  er 

Mow  with  the  help  of  A4.14  and  A4.12,  the  square  of  the  ratio, 
solid  particles  velocity  to  the  pressure  at  the  surface  can  be 
readily  computed  once  we  know  the  propagation  constants  k',  k" 
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I  In  the  case  of  normally  incident  (8q  »  0)  plane  wave,  the  ratio 

|  can  be  expressed  in  a  compact  form 


2 

..  2  2r 

V 

S 

(1-0)  w  |  1 

‘  i- 

s 

Lk7 

k"J 

P 

2.  2 

as  r  =  -1  and  =  0 


(A4.15) 


This  result  differs  substantially  from  that  obtained  by  Flohr  and 
Cress  (A4 . 1 Iwherein  the  ground  has  been  considered  as  an  elastic 
continuum.  The  present  result  is  frequency  dependent  which  is 
qualitatively  consistent  with  their  measurements. 
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The  effect  of  inhomogeneity  in  the  ground  on  the  wave  propagation 
can  be  described  by  considering  it  as  consisting  of  several  layers 
of  differing  properties.  One  such  simple  model  is  to  consider 
a  finite  depth  of  poro-elastic  layer  supported  by  a  rigid  semi- 
infinite  medium. 


In  this  section,  we  derive  an  analogous  expression  to  that  of 
Eq  A4.15  for  normally- incident  plane  waves.  When  such  a  plane 
wave  is  incident  upon  the  layered  medium,  inside  the  poro  elastic 
layer  of  depth  say,  d,  there  will  in  general  be  twelve  possible 
waves.  As  before,  they  can  be  represented  as  follows: 


In  the  solid  phase 

♦.  *  ”  a  '  exp  (-  iuit)  {e+  ik  n  +  T'e  ik  x  } 
IS 

♦  j"  «  as"  exp  (-  iuit)  { e+  ik"x  +T"e"  ilC^X  } 
ijfj  «  a#r  exp  (-  iuit)  te+  **  *  e  } 


(A4.16) 


In  the  fluid  phase 

♦2  pi 

*2"  *  T"  P2" 

*2  •TrP1 

The  unknown  quantities  t*  ,  T"  and  Tr  denote  the  relative  amplitudes 
of  backward  travelling  waves  to  the  waves  generated  at  the  surface. 


Now  at  the  rigid  backing  (at  x  ■  d)  the  tangential  stress  (o^)  will  be 
sero,  which  in  turn  from  stress-strain  relations  A4.7  implies  that 
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e^  =  0.  Furthermore  as  the  backing  is  hard  the  normal  particle 
velocity  is  zero  at  the  backing  i.e.. 


v  =  v  »  0  at  x  *  d 

s  t 


Also  at  the  surface  n  =  0,  the  tangential  stress,  e^  «  0. 


Using  the  above  conditions,  after  lengthy  algebra,  we  obtain  the 
following  for  the  square  of  the  ratio  of  solid  particle  velocity 


to  pressure  at  the  surface 


)2  =  (1  -  0) 2  u2 


[r(!  I  !*)'*•(!  1^)] 


2  2 
Kf“  <t*  -  t«r 


(A4.17) 


a  »  ik'd  and  6  «  lk"d 
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APPENDIX  5 

MOTION  OF  THE  INTERFACE.  THE  GEOMETRIC  APPROXIMATION 

Using  equations  <6.18)  and  the  Integral  expressions  for  the  potentials, 
and  after  performing  the  appropriate  differentiations  with  respect 
to  z,  we  obtain  the  components  of  displacement  at  the  interface  where 

Z  *  0: 


a_ 

3r 


-vizo 


7-  [a3+  aiaA2  -  V..A.J  - - Jo  (Kr)KdX 

Aq  U  -*  \)e 


Vi 


.  -ViZo 

[vsAi  +  a,,v2A2  ♦  K*Ai,]  ^-J0(Kr)KdK 


These  integrals  are  now  in  a  form  familiar  in  many  analyses  of  reflected 
spherical  waves  in  which  the  reflected  field  is  expressed  as  a 
superposition  of  plane  waves  multiplied  by  the  appropriate  reflection 
coefficient  (Brekhovskikh  [6.5]).  In  these  the  contribution  from 
the  saddle  point  (at  K  •  k}  sinB ,  where  8  is  the  angle  of  incidence) 
is  obtained  as  an  asymptotic  expansion,  the  first  term  of  which  gives 
the  geometric  approximation.  Ignoring  the  contributions  from  the 
branch  cuts  and  poles,  w«  have 


(A5.1) 


(A5.2) 


Ur  '  3r 


ikjR 


—  (A,  ♦  a,.  A,  -  (k*  sin*  8  -  )  *Afc»— -  U  ♦  OCR*1)) 

Aii  *  *  R 


L  a, 


— l-|~  Aj  +  a, ,  A2  -  <k\sin’  8  -  kl)JA»j  - 


ikiR 


(1  +  0(R-1)) 


(A5.3) 


S 


3: 

\  | 

si! 

B 

0 

I 

I 

I 


and 


u2  _  QkisinJ$  -  kJ)^A3  +  a,,  Casin’ 6  -  kJ)^A2 

ikiR 

+  k’jsin'e  Aw]  — R—  (1  +  0(R  )) 

where  R1  =  rJ  +  z\ ,  sin8  =  ~,  and  the  determinants  are  evaluated 
for  K  =  kjsinS. 

He  do  not  expect  the  geometric  approximation  to  be  valid  at  grazing 
incidence  because  the  determinants  A,  A,  Aw  all  contain  a  factor 

I  • 

vi ,  and  thus  approach  zero  as  t)  —>w/2.  To  see  this  we  recall  that 
they  were  formed  by  replacing  columns  2,  3,  or  4  of  the  matrix 
in  equation  (6.47)  by  the  vector  on  the  right-hand  side.  This  new 
column  may  then  be  added  to  column  1  without  affecting  the  value  of 
the  determinant ,  and  this  results  in  a  single  non-zero  element 
2v!  in  the  first  column,  which  must  be  a  factor  of  the  determinant. 

The  evaluation  of  the  higher  order  correction  terms  in  the  saddle  - 
point  approximation  would  be  over  complicated  because  of  the  need 
for  1st  and  2nd  derivatives  of  the  various  determinants.  However, 
we  later  use  a  simplified  method  which  enables  us  to  obtain  the 
correction  term  at  grazing  incidence. 


(A5.4) 
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A5 .  i  The  branch  cut  integrals 


.1 


if  we  treat  each  branch  cut  integral  separately,  then  our 
problem  in  general  is  to  evaluate  an  integral  of  the  form 
<* 

I  =  F(K)  Jg  (Kr)KdK 

<  O 

where  F(K)  has  a  single  branch  point  at,  say,  K  «  k_..  If  we 
write 

> 

F(K)e"Vj2J  (Kr)KdK 
o 

where  Vj?  ■  Rev  >o,  then  the  integration  around  the 

branch  cut  becomes  implicit. 


lim 

z'-*o 


We  now  expand  F(K)  as  a  Taylor  series  in  V^,  so  that 
F(K)  -  FdCj)  ♦  §  V3  ♦  i  ^  Vj  ♦ 

we  now  require  the  integrals 


lim 
z  '*o 


J0(Kr)KdX 


2 

where  R  ' 


♦  x 


•  0 


lim 
z  -ro 


V 


-v 


Jc(Kr)KdK 


lim 

z’-*o 


ik.r 


-I-  (IK)  -  ^  . 


i 
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r‘ 
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where  we  have  used  the  integral  representation  for  a  source  of 
spherical  waves,  as  in  section  6.8.  These  are  sufficient  to 
enable  us  to  evaluate  the  integral  to  terms  in  _1.  /  end  we  now  have 


I  -  /  dF 


iXie^S' 


+  0(r  ) 


Inspection  of  the  matrix  in  equation  (6.47)  shows  that  for  all 
terms  in  equations  (A5.1),  (A5.2)  we  can  write  in  general  form 


OO+VVh  (K) 


g2  (K)  +  Ujh,  (K) 


giving 


h  (k, )  g,  (ki  )-h,  (k-t )  g  (k-t  ) 


tg,  ocj  >  ia 


The  analysis  given  here  has  not  been  rigorous,  and  certainly  one 
questionable  item  is  the  convergence  of  the  Taylor  series  for  F(K) 
Without  giving  a  detailed  discussion,  this  will  be  in  doubt  when 
g2  (k._j  )  is  small,  implying  a  pole  near  the  branch  point.  To  find 
an  upper  limit  for  the  branch  cut  integral  we  ignore  g2 ,  leaving 

F(K)  «  g  (K)  ^  h  (K) 

_1 _  +  _J _ 

Vjh,  (K)  h,  (K) 


The  second  term  gives  no  contribution,  and  we  now  have 


I  ■  lim 


gi(K)  e-^Z' 

h,  (K)  V: 

O  J 


j„ (Kr)KdK 


Using  the  methods  of  section  6.11,  this  gives 


I  -  a  Otj  )  e  l  -* 

2,^42  ® — 3—  +  0(r  ) 
h,  Ot j  )  r 
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We  therefore  expect  equation  A5.5  to  be  valid  only  if 


<  «  ik  .c 

(kj)  e  j 


M*-)) 


» 


iki  (h1(jc3)g«(;y';h»(ltj)g1{lcj)) 

[g,  (*.,)  1 


lk  .r 


rg,(kj)  [g,  (kj)  1* _  >>  j  (A5.9) 

k  h,  (k^fh^k^g^k^)  -  VVW1 

This  quantity  may  be  regarded  as  related  to  the  numerical  distance 
established  more  precisely  for  the  rigid  porous  medium. 


A5-2  Grazing  incidence 

Although  the  evaluation  of  the  correction  term  in  the  saddle-point 
integral  is  impractical  for  an  elevated  source,  when  z.q  =■  0  the 
integral  may  be  handled  in  exactly  the  same  way  as  any  other  branch  cut 
integral.  This  gives  us  the  correction  term  for  grazing  incidence,  and 
this  may  be  a  useful  approximation  for  all  cases  where  the  geometric 
term  is  not  dominant. 


As  a  demonstration,  we  again  consider  the  field  at  the  surface  of 
rigid  porous  half-space  at  grazing  incidence.  This  is  given  by 
a  _  2i 


krie  (na-l) 


eikir 

®— +  0  (r  ) 


The  upper  limit  is 


a  ,  „ik  r 
v  -  2  e _ 


and  ($.64)  becomes 
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|kirM,(n*  -  1) 

These  results  are  the  limiting  case  of  the  more  elaborate  saddle-point 
integration. 


»  1 


A5.3  The  kj  branch  cut  integrals 

To  keep  the  expressions  for  the  integrals  as  compact  as  possible,  we 
define  V  as  the  cofactor  of  the  matrix  element  in  row  j  and  column  k 

1  jK 

in  the  expansion  of  the  determinant  A^.  He  then  find  expressions  for 
g^,  etc.  in  the  general  equation  (A5.6) 

When  dealing  with  the  ki  branch  cut  integral,  we  must  have  zq  ~  0.  In 
equation  (A5 . 1 )  we  obtain 

gi  (kj)  *=  2(6»n(ki)+  aj2A»xi  (ki)  “  (kj1  -to,1 )  ^61,11) 


hi(kx)  -  0 

9,  ( ki J  -  Ao{kx)  h,  (k^  «  ^,u(ki) 
and  in  equation  (A5.2) 


gitk!)  =  2  (<ki‘-k  *)*6ni  +  aM(ki,-K.,)i6gn+  kj’Sm) 
hj  (kj)  -  0,  g,  (kj)  -Ao(kj),  h,(ki)  mS0ll<ki' 


Using  equation 


2  ,  e1Xir 

3r  1  r1 


(A5.5), 

ikLe 


and  noting  that 
ik2r 

-  +0(r"’) 


r* 


we  obtain  the  displacements  in  the  far-field 
u  t  •  2  k  i*6o^[6in+  a»26tix~  (ki*  -  k»*)^5»n] 


A 

0 


2 

r 


(AS. 10) 
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A5.4  The  k2  branch  cut  inteqrals 

For  this  and  the  remaining  branch  cut  Integrals  we  are  not  limited 
to  z,  =0,  but  otherwise  we  proceed  as  before.  For  u^  we  have 


e“VlZo 

9i  =  — — —  (A,  +  a, ,4,  -  Vi»A*l 
e“VlZ# 

hi  *  — - -  *  *saAa»2  ~  Vn(S*i2  +  a324m,i)] 


g,  »  Ao 


h,  »  5oi2  +  a,j4o<.2 


(A5.13) 


(A5 .14) 


(A5.15) 


(A5.16) 


and  for  u 

z 

-ViZ0 

9l  “  — —  tv,&1  +  Aa> 1  (A5.17) 

-vizo 

hi  » — - -  [a, i  Aj  +  v>^«it  +  +  +  ®f  f  Ai| *2 ) )  (A5.18) 

l 

g2  *  Ao  /  h2  *  Aoj2  +  a32Ao%2  (AS.  19)  t  (A5.20) 


where  it  is  assumed  that  all  expressions  are  evaluated  for  K  *  k, 
We  shall  not  carry  out  here  the  straightforward  substitution  into 
equations  (A5.7)  and  (A5.5). 


A5.7  Surface  waves 


Surface  wave  contributions  will  occur  when  poles  of  the  various  integrands  lie 
between  the  real  axis  and  the  contour  of  integration.  For  this  to  happen 
we  must  have  =  0. 

When  the  wavenumbers  of  the  different  wave  types  are  real,  a  true  surface 
wave  must  have  a  real  wavenumber  greater  than  these  ([6.3],  [6.4]). 

When  at  least  one  wavenumber  is  complex,  as  in  our  application,  the 
situation  is  less  clear.  In  general  it  will  be  necessary  to  search 
for  zeros  of  A0  by  some  numerical  means,  and  also  to  establish  criteria 
as  to  whether  they  lie  within  the  closed  integration  contour. 

The  general  form  of  the  required  integral  is 

OP 

I  =  F(k)  J0  (Kr)KdX 

which  may  be  rewritten  as 
00 

I  =  1  F(K)hJ  (Kr)KdX  (A5.37) 

-flO 

In  the  simplifying  case  of  grasing  incidence  the  required  Integration 

contours  are  those  which  pass  around  the  vertical  branch  cuts  emanating 

from  the  branch  points  at  k,f  k_,  k, ,  k.. 

1  2  J  4 


A  typical  situation  is  shown  in  the  diagram,  where  poles  in  any  of  the 
regions  I  -  V  will  contribute  to  the  integral.  The  signs  of  the 
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In  region  III 
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In  region  V 
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Here  we  have  assumed  that  the  pole  lies  below  the  lashed  line  which 

consists  of  the  contours  Re  v  >0  and  In  e  >0. 

2  2 


We  now  write  the  integrals  as 
I 


I,  ♦  I,  ♦  I.  ♦  I.  +  2*i  res  (K  ) 

1  2  J  4  p 


where  I , ,I„  etc.  are  the  branch  cut  integrals,  anc  rae(K  )  represents 
l  i  p 

the  residues  of  any  poles  lying  within  the  prescribed  regions. 


(A5.40) 


I 

where  vl«v;>'vvv/t'  A  Q«  a2'  A y  A/)  are  evaluated  for  K  =  Kp. 

The  major  difficulties  are  the  determination  of  ,  and  the  evaluation 

I 

of  A  ^  (Kp )  .  In  the  next  section  we  use  a  light  fluid  approximation 
which  simplifies  the  problem  considerably. 


A5.8  The  light  fluid  limit. 

A  limiting  case  of  sane  Interest  occurs  when  the  density  of  the  fluid 
is  much  less  than  that  of  the  solid  in  the  lower  medium,  and  we  expect 
this  to  be  applicable  when  modelling  ground  surfaces.  In  this  case 
we  have  p,  p^  <<  p ,  and  also  recall  from  the  theory  of  rigid 

porous  media  that  Pj^k  c,  at  low  frequencies.  Since  we  expect  the 
propagation  constants  of  the  various  wave  types  to  be  roughly  comparable , 
then  the  compressibility  of  the  fluid  must  also  be  much  greater  than 


that  of  the  solid 


J 
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If  we  examine  the  effect  of  this  on  the  fast  and  slow  waves,  then  the 
approximations 


are  valid  if 


M 

H 


<< 


1. 


2 


It  should  be  remarked  here  that  a  is  always  less  than  unit  y  in  Biot's 
theory . 


For  the  shear  wave 


k  2  «  S-e  if 

4  g 


«  1 


PP. 


From  equation  6 . 36 
2  P 


32  P 


k'2 

p  v—  -  c0/<  —  -  1) 


-  a)/(l  -  _3 _  ) 

<2 


and  if  k  ^  <  k  ~  «  1  implies  a  <<  1 
J  2  H  32 


From  equation  6.37 


3  °i  3 

*23  “  -  ^)/d 

v  *2  2  v  '2 


(A5.41) 


(A5.42) 


(A5.43) 


These  expressions  sijq>lify  further  if  O^/Pj  K<  °>  and  *3  <<  kj  * 


Since  the  fast  wave  will  invariably  have  the  highest  velocity  of  all 
wave  types,  this  latter  condition  is  effectively  one  of  local 
reaction.  Under  these  circumstances  we  have 


a 


32 


-an 

H 


oM  p 

*23  =  H  p 


In  the  light  fluid  limit  we  can  also  show  that 


(A5.44) 


(AS.  45) 


A  -  A  A.  (A5.46) 

o  a  b 


where 


A^  =  u2  (P2V1  +  P3V  2^  (*5.47) 

Afa  -  (2pK2  -  kj^)  <>*2  -  K2) 

-  2|iv3V4K2  (A5.48) 

It  can  be  seen  that  equation (A5. 4 7)  is  the  denominator  found  for  the 
rigid  porous  case,  and  (AS.  4  8)  is  the  Rayleigh  denostlnator  for  a  semi- 
in  finite  elastic  half-space.  This  is  intuitively  appealing  since  it 
means  that  in  the  light  fluid  limit  wave  propagation  in  the  upper 
medium  is  identical  to  that  for  the  rigid  porous  half-space,  while  the 
Rayleigh  surface  wave  la  virtually  unaffected  by  the  presence  of 
the  fluid. 

Air-coupled  surface  waves  of  ths  Rayleigh  type  have  been  found 


experimentally  [5.2],  so  we  suspect  that  poles  corresponding  to  the 
zeros  of(A5.48)  will  be  of  importance.  In  the  surface  wave  expressions 
given  in  section(A5 . 7)  we  now  have 


A5.9  Practical  Calculations 


The  calculation  of  the  plane  wave  and  branch  cut  contributions  present 

no  difficulty,  but  to  obtain  any  surface  wave  contribution  we  need  to 

find  a  zero  of  A  and  its  derivative  at  this  point  (see  section  AS. 7) . 
o 

To  do  this  in  practice  we  have  used  a  numerical  technique  known  as 

Muller's  method  (described  in  reference  A5.5)  which  finds  a  root  of  an 

analytic  complex  function  by  using  quadratic  interpolation.  Three 

initial  guesses  are  used  to  obtain  the  first  quadratic,  and  the  nearest 

zero  of  this  provides  the  next  approximation  to  the  unknown  root,  and 

the  iteration  is  continued  until  the  required  accuracy  is  achieved. 

This  root  gives  us  the  surface  wave  number  Kp.  consequently  it  is  a 

simple  matter  to  find  the  derivative  A*  (K  )  by  calculating  the  gradient 

o  p 

of  the  last  fitted  quadratic. 

From  the  considerations  of  the  previous  section  we  expect  a  surface 

wave  of  the  Rayleigh  type  which  will  have  a  wave  number  approximately 

1.2  times  that  of  the  shear  wave ,  k^.  Using  three  initial  guesses 

1.15k4,  1.2k4  and  1.25k4  resulted  in  rapid  convergence  for  all  the 

exaaples  tried,  and  there  was  no  indication  of  any  problem  arising  from 

the  non -ana lytic  nature  of  A  (due  to  the  branch  cuts) 

o 

AS. 10  Material  parameters  for  soil. 

Sound  velocity  (wA^)  3.4  x  10*  cm  s-1 

Density  of  air  (p^)  1.2  x  10  ^  g  cm  * 

p  -  wave  velocity  5.0  x  10*  cm  s"1 


J 


l 


im  <kj)/  Re(kj) 

1.25  x  10-2 

S  -  wave  velocity  (w/ReOc^) ) 

3.2  x  104  cm  ; 

im(k4)/Re(k4> 

1.25  x  10-2 
-3 

Density  of  porous  medium  (p) 

1.3  g  cm 

Tortuosity 

2.5 

Shape  factor 

0.75 

Flow  resistivity 

250  cgs  rayls 

Porosity 

0.4 

a 

0.5 

The  imaginary  parts  of  the  compressions!  and  shear  wave  numbers  are 
consistent  with  the  results  found  by  Prange  [AS. 6]  for  attenuation  in 
soils. 


The  tortuosity,  shape  factor,  flow  resistivity,  and  porosity  are  used 
to  determine  the  equivalent  rigid-porous  wave  number  kj  end  complex 
density  pj  (see  chapter  5).  The  value  chosen  for  the  parameter  a  is 
arbitrary,  apart  from  the  constraint  in  Blot's  theory  that  it  must  lie 
between  the  porosity  and  unity.  However,  we  did  repeat  some  of  the 
calculations  for  a  »  0.9,  end  found  that  the  changes  to  the  computed 
surface  velocity  were  not  more  than  2dB,  and  in  most  cases  considerably 
lass. 

A5.ll  Results 

Figures  A5.1  -  A5.S  show  the  contributions  of  the  various  wave  types  to 
the  vertical  component  of  surface  velocity  for  a  source  height  of  1  s 
and  a  horisontal  separation  of  100  m.  These  are  the  plane  wave 


I 


J 


approximation ,  the  n*u  and  far- fie  Id  approx  imat  ions  to  the  k^ ,  kj  and 
k^  branch  cut  integrals,  and  the  Rayleigh-type  surface  wave.  The  k^ 
branch  cut  integral  was  always  negligible  because  of  rapid  attenuation 
of  the  slow  wave. 


Figures  AS. 6  -  A5.8  show  the  dominant  contributions  for  source  heights 
of  zero,  1  d  and  SO  ■. 

In  all  these  cases  a  point  source  of  spherical  waves  is  used,  such  that 
the  pressure  amplitude  of  the  incident  wave  is 

ik<  R 
P  R 

The  results  show  that  with  the  source  on  the  surface  the  main 
contributions  are  those  of  the  Rayleigh-type  wave  and  the  k^  branch 
cut  integral.  The  Rayleigh  wave  dominates  below  200  Ha,  but  is 
attenuated  at  higher  frequencies.  Elevation  of  the  source  to  1  a 
gives  a  significant  plane  wave  contribution,  which  becomes  increasingly 
dominant  at  high  frequencies.  The  plane  wave  ten  dominates  at  all 
frequencies  when  the  source  height  is  increased  to  50  m.  There  are 
contributions  from  the  k ^  and  k^  branch  cuts,  but  these  are  always 
a  small  part  of  the  total  motion. 
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J  Figure  AS . 1  Plane  Nave  Contribution  for  a  Source  Height  of  1 B 
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Figure  A5.2  Contribution  from  k  Branch  Cut  for  a  Source  Height  of  In 

Near  field - L - 

Far  field  —  .  -  .  - 


188 

Figure  A5.3  Contributions  froa  k,  Branch  Cut  for  a  Source  Height  of  la 
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Figure  A5.4  Contribution  fro*  k^  Branch  Cut  for  a  Source  Height  of  1  a 
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Figure  AS. 5  Contribution  from  Rayleigh-type  Wave  for  a  Source  Height  of  lm 
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I  j  figure  A5.6  Contributions  from  ^  Cut  (Near-  and  Far-Field)  and  Rayleigh-Type  wave  for 
*  Source  and  Receiver  Both  on  Surface 
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Figure  AS. 7  Contributions  fro*  Plane  Nave,  k  Brandt  Cut  (Near-  and  Par-Field)  and 
Rayleigh-type  Nave  for  Source  Height  of  la 
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Figure  A5.8  Dominant  Plane-Have  Contribution  for  a  Source  Height  of  50m 
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